Geometry and dynamics of a tidally deformed black hole 
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The metric of a nonrotating black hole deformed by a tidal interaction is calculated and expressed 
as an expansion in the strength of the tidal coupling. The expansion parameter is the inverse length 
scale TZr^ , where TZ is the radius of curvature of the external spacetime in which the black hole 
moves. The expansion begins at order TZT^ , and it is carried out through order 7?.~*. The metric 
is parameterized by a number of tidal multipole moments, which specify the black hole's tidal 
environment. The tidal moments are freely-specifiable functions of time that are related to the 
Weyl tensor of the external spacetime. At order TZ"^ the metric involves the tidal quadrupole 
moments £ab and Bab- At order TZ~'^ it involves the time derivative of the quadrupole moments 
and the tidal octupole moments Eabc and Babe- At order 7?."'* the metric involves the second time 
derivative of the quadrupole moments, the first time derivative of the octupole moments, the tidal 
hexadecapole moments Eabcd and Babcd, and bilinear combinations of the quadrupole moments. The 
metric is presented in a light-cone coordinate system that possesses a clear geometrical meaning: 
The advanced-time coordinate v is constant on past light cones that converge toward the black hole; 
the angles 9 and (j> are constant on the null generators of each light cone; and the radial coordinate 
r is an affine parameter on each generator, which decreases as the light cones converge toward the 
black hole. The coordinates are well-behaved on the black-hole horizon, and they are adjusted 
so that the coordinate description of the horizon is the same as in the Schwarzschild geometry: 
r = 2M + 0{H~^). At the order of accuracy maintained in this work, the horizon is a stationary 
null hypersurface foliated by apparent horizons; it is an isolated horizon in the sense of Ashtekar 
and Krishnan. As an application of our results we examine the induced geometry and dynamics of 
the horizon, and calculate the rate at which the black-hole surface area increases as a result of the 
tidal interaction. 



I. INTRODUCTION AND OVERVIEW 
This work and its context 

Our main goal in this work is to calculate the gravita- 
tional field of a nonrotating black hole that is deformed 
by a tidal interaction with external bodies. We assume 
that the tidal interaction is weak, and that it changes 
slowly; we make no other assumptions, and describe the 
tidal environment in the most general terms compatible 
with the main assumptions. This work is a continuation 
of a line of inquiry that was initiated by Manasse [1^ in 
the early nineteen sixties; our implementation is much 
more general, and much more accurate, than Manasse 's 
original work. 

We first introduce the two length scales that are rel- 
evant to this problem. (Throughout the paper we use 
relativistic units and set G = c = 1.) The first is M, the 
mass of the black hole. The second is TZ, the radius of 
curvature of the external spacetime associated with the 
external bodies, evaluated at the black hole's position. 
We assume that the scales are widely separated, so that 

Af<7^. (1.1) 

This condition ensures that the tidal interaction is weak, 
and allows us to speak meaningfully of a black hole mov- 
ing in an external spacetime; when M is comparable to 



TZ, no clear distinction can be made between the "black 
hole" and the "external spacetime." We refer to the ap- 
proximation scheme based on the condition M/TZ ^ 1 as 
the small-tide approximation. 

As a concrete example we may consider a situation in 
which the black hole is a member of a binary system. Let 
M' denote the mass of the external body, h the separation 
between the companions, and V - ^{M + M')/b the 
orbital velocity. The radius of curvature is then TZ 
^b^/{M + M'), and 



TZ M + M' ■ ' ■ ' 

We demand that this be a small quantity. There are two 
particular ways to achieve this. In the small-hole ap- 
proximation the black-hole mass is assumed to be much 
smaller than the external mass, so that M/{M-\-M') <C 1; 
then M/TZ is small irrespective of the size of V , and 
the binary system can be strongly relativistic. In the 
weak-field approximation it is V that is assumed to be 
small, while the mass ratio is left unconstrained; here 
the two companions can have comparable masses, or the 
black hole can be much larger than its companion, but 
the mutual gravity between the bodies must be weak. 
The small-hole and weak-field approximations are par- 
ticular instances of the more general requirement that 
M/TZ <^ 1; they are both incorporated within our small- 
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tide approximation. 



Tidal quadrupole moments 



The effects of a tidal field on the structure of spacetime 
around a black hole were first investigated by Manasse 
[1], in the specific context of the small-hole approxima- 
tion. Using techniques similar to those exploited in this 
paper, Manasse calculated the metric of a small black 
hole that falls radially toward a much larger black hole. 
Each black hole was taken to be nonrotating, and the 
small hole was taken to move on a geodesic of the (un- 
perturbed) Schwarzschild spacetime of the large hole. 
The case of circular motion around a large Schwarzschild 
black hole was treated much later by Poisson [2], and 
Comeau and Poisson [5] examined the case of circular 
motion around a Kerr black hole. 



The methods employed by Manasse could be applied 
beyond the small-hole approximation. Alvi [HIS] realized 
that they could be seamlessly extended to the more gen- 
eral context of the small-tide approximation of Eq. ( |1.1 1 . 
Alvi exploited this insight to calculate the tidal fields 
acting on a black hole in a post-Newtonian binary sys- 
tem. In Alvi's work, the two bodies have comparable 
masses and the black hole has a significant influence on 
the geometry of the external spacetime. Alvi calculated 
the tidal fields to leading (Newtonian) order in the post- 
Newtonian approximation to general relativity, and spe- 
cialized the orbital motion to circular orbits. His work 
was later generalized to first post-Newtonian order, and 
to generic orbits, by Taylor and Poisson J6j. 



In the work of Detweiler (T) [5] and Poisson [S] reviewed 
previously, and in the work presented here, the tidal en- 
vironment is described in the most general terms com- 
patible with the Einstein field equations. The metric 
is parameterized by freely-specifiable functions of time 
that serve the specific purpose of specifying the tidal 
environment; these are packaged in symmetric tracefree 
(STF) tensors Sa^az - ai and Kaioa -ai llSj, which we re- 
fer to as the tidal multipole moments of the black-hole 
spacetime. (The tidal moments can be thought of as 
Cartesian tensors; they are symmetric and tracefree in 
all pairs of indices.) At order (r/TZ)'^ the metric in- 
volves the quadrupole moments Sab (5 functions) and 
Bab (5 functions). At order (r/TZ)^ the metric involves 
the time derivative of the quadrupole moments, as well 
as the octupole moments Sabc (7 functions) and Babe (7 
functions). And at order (r/TZ)'^ the metric involves the 
second time derivative of the quadrupole moments, the 
first time derivative of the octupole moments, the hex- 
adecapole moments Sabcd (9 functions) and Babcd (9 func- 
tions), and bilinear combinations of the quadrupole mo- 
ments. This generality is a substantial virtue of our work, 
which unifies and extends earlier works [H IH] [THl 120] that 
examined special cases. 



Alvi's work motivated an effort to improve our under- 
standing of the tidal interaction of black holes by con- 
structing the black-hole metric to high order in the cou- 
pling strength M/TZ. More precisely, the metric is cal- 
culated in the black hole's local neighborhood, and ex- 
pressed as an expansion in powers of r/7?. <C 1, where r 
is the distance to the black hole; the metric is valid to all 
orders in M/r. This program was initiated by Detweiler 
[7JIH], who calculated the metric through order (r/TZ)^, 
the lowest order at which tidal effects appear. It was pur- 
sued by Poisson [S], who calculated the metric through 
order (r/TZ)^. Here we calculate the metric through or- 
der (r/TZ)^, and present the most accurate version ever 
produced of the metric of a tidally deformed black hole. 

Other authors have also contributed to this effort. 
Frolov and his collaborators [lOl [11] examined the in- 
ternal geometry of a tidally deformed black hole, and 
Damour and Lecian [12; characterized the tidal deforma- 
tion in terms of a polarizability and "shape Love num- 
bers;" these works were restricted to axisymmetric tidal 
fields, a restriction that is not made in this paper. In 
a genuine tour de force, Yunes and Gonzalez |13j calcu- 
lated the tidal deformation of a rapidly rotating black 
hole to order {r /TZ)'^ . The tidal deformation of neutron 
stars (and other types of compact bodies) has also been 
the subject of recent investigations [21 [TSJ [111 117] ■ 



The metric of a tidally deformed black hole is obtained 
by integrating the vacuum field equations in the local 
neighborhood of the black hole. The field equations leave 
the tidal moments undetermined, and the metric is pre- 
sented as a functional of these arbitrary moments. In ap- 
plications of this framework the tidal environment must 
be specified, and this is achieved by making appropri- 
ate choices for the tidal moments. In practice this typi- 
cally requires matching the local metric to a global met- 
ric that includes the black hole and the external bodies 
that are responsible for the tidal interaction. For exam- 
ple, Taylor and Poisson [Bj and Johnson-McDaniel et al. 
[21] carried out such a matching in the context of the 
slow-motion approximation, in which the mutual gravity 
between the black hole and the external bodies is weak. 
The global metric was expressed as a post-Newtonian 
expansion, within which the black hole can justifiably 
be represented as a point particle. These authors de- 
termined the quadrupole moments £ab and Bab to first 
post-Newtonian order, and their work could easily be ex- 
tended to obtain the higher multipole moments that also 
appear in the black- hole metric. In general the nonlinear- 
ities of the field equations imply that tidal fields depend 
on the black-hole mass M. In the small-hole approxi- 
mation, however, the black hole can be treated as a test 
body, and this dependence disappears; in this case the 
determination of the tidal moments is simplified [T] [51 [3] . 



3 



Light-cone coordinates 

The choice of coordinates is often critical in the con- 
struction of a metric and the exploration of its properties. 
This is all the more true in the case of black-hole space- 
times, which require coordinates that are well-behaved 
on the horizon. We have given a lot of attention to 
the selection of coordinates, and have chosen to work 
with a system {v,r,d,(j>) that is specifically tailored to 
describe the geometry of past light cones. We refer to 
these as light- cone coordinates. The coordinates have 
a clear geometrical meaning: The advanced-time coor- 
dinate V is constant on past light cones that converge 
toward the black hole; the angles 9 and cj) (which we col- 
lectively denote 9^, with the upper-case Latin index A 
running from 2 to 3) are constant on the null genera- 
tors of each light cone; and the radial coordinate r is an 
affine parameter on each generator, which decreases as 
the light cones converge toward the black hole. Through 
order (r/TZ)^ the radial coordinate doubles as an areal 
radius, in the sense that the area of each two-surface 
{v,r) = constant is equal to 47rr^[l -I- 0{r'^ /TZ'^)]; this 
property is lost at order (r/TZ)'^. In addition, the radial 
coordinate is tuned so that the coordinate description of 
the black-hole horizon is the same as in the Schwarzschild 
spacetime: r = 2M[1 + 0{M^/n^)]. The light-cone co- 
ordinates are well-behaved across the horizon. 

The choice of coordinates is inspired from the work 
of Bondi et al. ^\ and Sachs [55], in which light- 
cone coordinates were utilized to construct the metric 
of an asymptotically-flat spacetime. (Here the coordi- 
nates were based on future light cones that expand to- 
ward future null infinity.) It is inspired also by the work 
of Ellis and collaborators on observational cosmology 
[211251 1571 [251 HH] , in which the metric of an expand- 
ing universe is constructed from observations made by 
a typical cosmological observer. Light-cone coordinates 
were also introduced by Synge in his classic textbook on 
general relativity [3D] ; we followed his methods closely in 
this work and its precursors [21 1^ [511 152] . 

In addition to the quasi-spherical system of light-cone 
coordinates {v,r,9^), we find it convenient to introduce 
also a quasi-Cartesian variant {v,x°'). The spatial coor- 
dinates x° (with the lower-case Latin index a running 
from 1 to 3) are constructed in the usual way from the 
quasi-spherical coordinates {r,9^); we have the relations 
X = r sin 9 cos (j), y = r sin 9 sin 0, and z = r cos 9, which 
we collectively denote x°- = rri°(0'^). 



preserves the geometrical meaning of the background co- 
ordinates. To construct the perturbation we rely on the 
covariant and gauge-invariant formalism of Martel and 
Poisson [33], and on the formulation of the light-cone 
gauge by Preston and Poisson [31]. At orders {r/TZY 
and {r/TZY the perturbation satisfies the vacuum field 
equations linearized with respect to the Schwarzschild 
solution. At order (r/TZ)'^ the perturbation satisfies non- 
linear field equations. 

In the quasi-Cartesian coordinates {VjX"") the build- 
ing blocks of the metric are the tidal potentials intro- 
duced in Tables |!j|l!| [m] |!vl |V| and |Vl]of Sec. |ll| these 
are generated by the tidal multipole moments introduced 
previously. The black-hole metric appears in Eqs. (3.5 1 
of Sec. Ill C| and it is expressed in terms of the radial 
functions listed in Tables jXIVj and |XV[ In the quasi- 
spherical coordinates (w,r, 6*"^) the tidal potentials are 
listed in Tables \Vm\ |IXJ |X] [XT] |XlT] and [XIIIl of Sec. [IT] 
and expressed as expansions in spherical-harmonic func- 
tions (see Table Vlll. The black-hole metric appears in 



Eqs. (3.7) of Sec. IIIC and it involves the same set of 
radial functions. 

An important property of the black-hole metric is the 
fact, mentioned previously, that in our light-cone coordi- 
nates, the position of the horizon is given by 



rho 



= 2M[1 + 0(MV7^^)], 



(1.3) 



the same relation as in the unperturbed Schwarzschild 
spacetime. This information allows us to investigate the 
nature and dynamics of the horizon's geometry, which 
is determined by the horizon's induced metric as well as 
the expansion and shear of the horizon's generators. An 
important outcome of this investigation is the statement 
that 



(1.4) 



This represents the rate at which the congruence of null 
generators expand, and we find that within the degree 
of accuracy maintained in this work, the generators are 
stationary. This implies that the horizon of a tidally 
deformed black hole is foliated by apparent horizons. The 
black-hole horizon is an isolated horizon in the sense of 
Ashtekar and Krishnan [351 IMl ISZ] • 



Tidal heating 



Black-hole metric 

The metric of a tidally deformed black hole is obtained 
by constructing a perturbation of the Schwarzschild so- 
lution, which describes a nonrotating black hole in com- 
plete isolation. The Schwarzschild metric is presented 
in the Eddington-Finkelstein coordinates {v,r,9"^), and 
the perturbation is presented in a light-cone gauge that 



A dynamical consequence of the tidal interaction is the 
fact that the black hole grows in size at a rate described 

by 



+ o(MV7^9). 



(1.5) 
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Here A is the surface area of a cross-section v = constant 
of the black- hole horizon, and kq := (4M)~^ is the sur- 
face gravity of the unperturbed horizon. An overdot indi- 
cates differentiation with respect to v, and the right-hand 
side of the equation involves the tidal moments Sab, £abc, 
Bab, and Babe- We refer to the growth of area that results 
from the tidal interaction as the tidal heating of the black 
hole by the external bodies. This choice of terminology 
deserves an explanation. 

The changes in black-hole parameters that result from 
time-dependent, external processes have been the subject 
of investigation by many researchers, starting from the 
pioneering work of Teukolsky and Press [3S] . It is useful 
to classify the physical processes that alter the configu- 
ration of a black hole as fast processes on the one hand, 
and slow processes on the other. Fast processes occur on 
a time scale that is shorter than or comparable to the 
black-hole mass M, and these produce (electromagnetic 
and/or gravitational) radiation that is partially absorbed 
by the black hole; in this case the changes in black-hole 
parameters that result from the interaction are radiative 
changes. Slow processes, on the other hand, occur on 
a time scale that is long compared with the black-hole 
mass; while these continue to produce changes in the 
black-hole parameters, the phenomenon no longer pos- 
sesses a radiative character, and the phrase "tidal heat- 
ing" captures the physics better than the phrase "black- 
hole absorption." It is good to point out that the math- 
ematical formulation of the phenomenon by Teukolsky 
and Press [38^ (see also Ref. [39^) is valid both for fast 
and slow processes, and is insensitive to matters of inter- 
pretation. 

The tidal heating of black holes was recently investi- 
gated by Alvi gO], Poisson Taylor and Poisson [5], 
and Comeau and Poisson [3], building on earlier work by 
Poisson and Sasaki [41] and Tagoshi, Mano, and Taka- 
sugi |42|. The notion of tidal work, tidal torque, and 
tidal heating was put on a firm relativistic footing by 
Purdue [43 , Favata [H], and Booth and Creighton |15] . 
In a recent work [46], Poisson compared the equations 
that describe the rate of change of the black-hole mass, 
angular momentum, and surface area that result from a 
tidal interaction with external bodies, with the equations 
that describe how tidal forces do work, torque, and pro- 
duce heat in a Newtonian, viscous body; the equations 
are strikingly similar, and the correspondence between 
the Newtonian-body and black-hole results is revealed to 
hold in near-quantitative detail. The tidal heating and 
torquing of black holes was incorporated in an effective 
theory of point particles by Goldberger and Rothstein 
[Hj and Porto [48 j. 

In favorable circumstances the tidal heating and 
torquing of a black hole can be relevant to astrophysi- 
cal sources of gravitational waves US] . In particular, it is 
likely to be significant in the generation of low-frequency 
waves that would be measured by a space-based detec- 
tor such as LISA [5^. For example, Martel [FT] showed 
that during a close encounter between a massive black 



hole and a compact body of a much smaller mass, up to 
approximately five percent of the lost orbital energy goes 
toward the tidal heating of the black hole; the rest is car- 
ried off to infinity by the gravitational waves. Hughes [52] 
calculated that when the massive black hole is rapidly ro- 
tating, the tidal heating slows down the inspiral of the 
orbiting body, and therefore increases the duration of 
the gravitational- wave signal. These conclusions are sup- 
ported by Hughes et al. [53] and Drasco and Hughes [53] . 

There are indications that the tidal heating and 
torquing of black holes may have been seen in accurate 
numerical simulations of the inspiral and merger of bi- 
nary black holes ^55] [55]. It is conceivable that as the 
precision of these simulations continue to improve in the 
future, the tidal heating will be exploited as a diagnos- 
tic of numerical accuracy; the surface area of each simu- 
lated black hole should be seen to grow in accordance to 
Eq. (1.5 1 instead of staying constant in time. 



Other applications 

The work presented here can serve as a foundation for 
the construction of initial data sets for the numerical sim- 
ulation of black-hole inspirals (see Ref. [57^ for a review 
of the current state of the art). Simulations carried out 
thus far have largely relied on initial data [581 ISS] that 
were adopted more for their flexibility and convenience 
than their astrophysical realism; these initial data tend 
to contain a large amount of spurious radiation and pro- 
duce unwanted eccentric orbital motion. A promising al- 
ternative strategy is to rely on a post-Newtonian metric 
to construct initial data sets that give a faithful (though 
approximate) description of two widely separated black 
holes. Early implementations [OQj |6TJ [62] of this idea 
treated the black holes as post-Newtonian point masses, 
and produced initial data that were unreliable close to 
each body. An improvement was proposed by Alvi [T, 'S], 
who recognized that the post-Newtonian metric should 
be replaced by a different representation of the gravita- 
tional field in the local neighborhood of each black hole; 
he therefore patched the post-Newtonian metric to the 
metric of a tidally deformed black hole in a buffer region 
M <^ r ^ TZ, in which each metric gives an acceptable 
representation of the true gravitational field. Alvi's con- 
struction was perfected by Yunes and his collaborators 
|63l 1641 165] . and the most mature implementation of this 
idea is contained in the recent work of Johnson-McDaniel 
et al. [H] . Because the black- hole metric presented in this 
paper is more accurate than the Detweiler metric [T] |8] 
used by Johnson-McDaniel et al., it could be involved in 
an improved version of their construction. 

Another interesting avenue of application for our met- 
ric would be to involve it in the dynamical-horizon for- 
malism of Ashtekar and Krishnan [351 IMl [22] , which pro- 
vides a purely local characterization of the structure and 
dynamics of black-hole horizons. Our work could con- 
tribute new insights to this effort by allowing the horizon 
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quantities to be expressed in terms of the tidal moments, 
which encode information about the external universe. 
In this way, our metric would be play the role of mes- 
senger between the outside world and the horizon. Ini- 
tial steps along those lines were taken by Kavanagh and 
Booth [66j . A particularly interesting question to inves- 
tigate is whether the dynamical-horizon notions of mass 
and current multipole moments |67j are compatible with 
the recent observation that the tidal Love numbers of a 
nonrotating black hole must all be zero [17]. 



Organization of the paper 



The remainder of the paper is divided into five sections 
and four appendices. We describe our results in Sees. [TTl 
|III| and |IV| and provide derivations of these results in 
Sees. |V] |VI| and the Appendices. 

We begin in Sec.|TT]by providing definitions for the tidal 
multipole moments Sab, £abc, £abcd, Bab, Babe, and Babcd- 
The tidal moments allow us to introduce the length scale 
TZ, and they give rise to tidal potentials that form the 
building blocks for the construction of the black-hole 
metric. 



In Sec. |III| we display our expressions for the metric of 
a tidally deformed black hole. We proceed in two steps. 
We first consider a smooth timelike geodesic 7 in a vac- 
uum region of an arbitrary spacetime, and we construct 
the metric of this spacetime in a neighborhood of the 
world line. We refer to this spacetime as the background 
spac etime, and its metric is displayed in Eqs. (3.3 1 and 
(3.4); the metric is presented in the light-cone coordi- 
nates {v,r,9^) and (w,a;°) introduced previously. We 
next insert a black hole of mass M into the background 
spacetime, place it on the world line 7, and re-calculate 
the metric. The metric of the black-hole spacetime is 
displayed in Eqs. (3.7 1 and (3.51. 



In Sec. |IV] we describe the consequences of the tidal 
deformation on the structure and dynamics of the black- 
hole horizo n. W e establish the statements of Eqs. (1.3 1, 
(fOl), and pTsl. 



metric of Eqs. (3.3 1 and (3.4 1. In Sec. VI we present a 



derivation of the black-hole metric of Eqs. ( 3.5 ) and ( 3.7 1 



In Appendix we describe how the tidal potentials 
can be decomposed in (scalar, vector, and tensor) spher- 
ical harmonics. In Appendix [B] we calculate the deter- 
minant of the induced metric on the black-hole horizon; 
the result is displayed (but not derived) in Sec. IV In 



Appendix [C] we provide calculational details relevant to 
the computation of the tidal heating in Sec. |IV| And in 
Appendix[D]we sketch a second, alternative derivation of 
the background metric of Eqs. (3.3 1 and (3.4 1. 



II. TIDAL MOMENTS AND POTENTIALS 

The work of Zhang [18] reveals that the metric of any 
vacuum spacetime can be constructed in the neighbor- 
hood of any geodesic world line and expressed in terms 
of two sets of tidal multipole moments. We begin our 
exploration of the geometry of a tidally deformed black 
hole in Sec. Ill Al with a formal definition of these mo- 
ments. We use them in Sec. |IIB| to specify the length 
and time scales associated with the tidal environment, 
and in Sec. |IIC] we describe their properties under parity 
transformations. We conclude in Sees. Ill D] and UTeI with 
the introduction of potentials that can be constructed 
from the tidal moments. 



A. Definition of tidal moments 

We consider a vacuum region of spacetime in a neigh- 
borhood of a smooth timelike geodesic 7. The world line 
is described by the parametric relations z"(t) in an ar- 
bitrary coordinate system x", and it is parameterized by 
proper time r. The velocity vector u" :— dz'^ /dr is tan- 
gent to the world line, and we construct a vectorial basis 
by adding to an orthonormal triad of vectors e"(r), 
which we assume to be orthogonal to and parallel- 
transported on the world line; the Latin index a labels 
the three members of the triad. 

We use the basis (u",e") to decompose tensors that 
are evaluated on the world line. For example. 



CabcO -^C^, .".^.T.M 



Cabcd ■ — Cafi-fSGaG^S^efi 



(2.1a) 
(2.1b) 
(2.1c) 



are the frame components of the Weyl tensor evaluated 
on 7; these are functions of proper time r. We shall also 
need the frame components of its first covariant deriva- 
tive, 



Cai 



,0&0|c 



abcQ\d 
abcd\e 



— ^ Oi^pv\^^a ^b c, 

— Cap-tii-se-ae^&Z^d^^ ^ 

— C ^ r P°' P^ p'^ P^ p'^ 



In Sec. [V] we provide a derivation of the background as well as 



cx fj, f3 '*y (5 
aQbO\cd '■— '^afi/3iy;jSSaU Sj^U 6^6^, 

CabcO\de '■— C af}^ fi-^Se^a^b^'c^^ ^d^e, 



c, 



ahcQ\de 
ahcd\ ef 



f^pip^.p>^p-^ 



al3-f5;eri'ia ^b ^c ^d'^e^f ) 



(2.2a) 
(2.2b) 
(2.2c) 

(2.3a) 
(2.3b) 
(2.3c) 



the frame components of its second covariant derivatives. 
We manipulate frame indices as if they were associated 
with Cartesian tensors; we lower them with the Kro- 
necker delta Sab, and we raise them with 6°^^. 

The symmetries of the Weyl tensor imply that it 
possesses 10 algebraically-independent components, and 
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these can be encoded in the two symmetric-tracefree 
(STF) tensors 



STF 



^ab (C'aObo) 
Bab ■= 2 (^apgC'^'bo) 



(2.4a) 
(2.4b) 



Here tabc is the permutation symbol, and the STF sign 
instructs us to symmetrize all free indices and remove all 
traces. We also use an angular-bracket notation to indi- 
cate the same operation: if Aabcd is an arbitrary tensor, 
then 



^(abcrf) (Aabcd) 



STF 



(2.5) 



In the case of Eq. (2.4 1 the STF operation is superfluous, 
because Caobo and iapqC^\o are already symmetric and 
tracefree. We refer to Sat and Bab as the tidal quadrupole 
moments associated with the world line 7. Each STF 
tensor contains 5 independent components, and these are 
functions of proper time r. We use Sat and Bab to denote 
the derivative of each moment with respect to r, and Sab 
and Bab to denote the second derivatives. The derivatives 
of the Weyl tensor in the direction of u" can be expressed 
directly in terms of these quantities. 

The symmetries of the Weyl tensor and the Bianchi 
identities imply that the spatial derivatives of the 
Weyl tensor — those listed in Eq. (2.2 1 — possess 24 
algebraically-independent components. These are en- 
coded in the STF tensors 



STF 



Sabc '■— (CaOf)0|c) 
Uabc •— g yi^apq'-J f,0|c/ ' 



(2.6a) 
(2.6b) 



and in Sab and Bab ■ We refer to Eabc and Babe as the tidal 
octupole moments associated with the world line 7. Each 
tensor contains 7 independent components, and these are 
functions of proper time r. We use Sabc and Babe to 
denote the derivative of each moment with respect to t. 
The derivatives of the spatially-differentiated Weyl tensor 
in the direction of can be expressed directly in terms 
of these quantities. 

We shall also require the secon d sp atial derivatives of 
the Weyl tensor, as listed in Eqs. (2.3 1. In this case there 



are 62 algebraically-independent components, and these 
are encoded in Sab, Bab, Sab, Bab, Sabe, Babe, and in the 
new STF tensors 

^abcd '■= ^{CaObO\cd) , (2.7a) 
Babed ^{eapqCZoledr'' ■ (2.7b) 

We refer to Sabcd and Babcd as the tidal hexadecapole mo- 
ments associated with the world line 7. Each tensor con- 
tains 9 independent components, and these are functions 
of proper time r. 

The decomposition of the Weyl tensor and its deriva- 



tives in terms of tidal moments is displayed in Eqs. (D15 I 



(D17), and (D19l of Appendix |D The numerical factors 
inserted in Eqs. (2.4 1, (2.6 1, and (2.7 1 are inherited from 



Zhang's choice of normalization |18j . 

The black-hole metric of Sec. |III| will be expressed in 
terms of the tidal moments Sab, Sabe, Sabcd, Bab, Babe, and 
Babed, which are treated as freely-specifiable functions 
of time. As we shall explain in Sec. |III| the relation 
between the tidal moments and the Weyl tensor will be 
more subtle than what was described here. 



B. Tidal scales 

The tidal moments allow us to specify length and time 
scales that characterize the tidal environment around the 
world line 7. We thus introduce the length scales TZ and 
C, the time scale T, and the velocity scale V. Our dis- 
cussion here follows Thorne and Hartle [68' and Zhang 

m- 

The length TZ is the local radius of curvature, and rep- 
resents the strength of the Weyl tensor evaluated on the 
world line. We express this as 



Sab 



1 



Bab 



7^2■ 



(2., 



The first equation indicates that a typical component 
of Sab will have a magnitude comparable to TZ~^. The 
second equation indicates that a typical component of 
Bab will differ from this by a factor of V, and this defines 
the velocity scale. 

The length C is the inhomogeneity scale, and it mea- 
sures the degree of spatial variation of the Weyl tensor. 
It is defined through the relations 



Sab 



1 



Bab 



V 



(2.9) 



involving the tidal octupole moments. We expect that 
the hexadecapole moments will be suppressed by an ad- 
ditional factor of C: 



Sa 



1 



bed 



Ba 



V 



bed 



(2.10) 



The time T is the scale associated with changes in the 
behavior of the Weyl tensor. This is defined by 



Sab 



Bab 



V 

TVT' 



(2.11) 



-1 



We also have Sab = 'R-'^T^'^ and Sabe ^ 7^"2£-l'- 
as well as Bab = VW^T-^ and Babe = VRr^C-^T- 

To illustrate the meaning of these tidal scales, let us 
consider as an example the world line of an observer mov- 
ing on a circular orbit of radius b around a body of mass 
M' . In this case the velocity scale is V ~ ^J~M^~lh, the 
radius of curvature is 7?. ^ yjlfi jM' , the inhomogeneity 
scale is £ ~ 6, and the time scale is T yJWjKT . In this 
example we have that T ^ TZ, V ^ -C/T, and C ^ VTZ. 
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TABLE I; Irreducible tidal potentials: type-f , even-parity. 
The superscripts q, o, and h stand for "quadrupole," "oc- 
tupole," and "hexadecapole," respectively. The tidal mo- 
ments £ab, Sabc, and £abcd are the STF tensors defined in 
the text. The extra factors of 2 in the hexadecapole poten- 
tials are inserted to respect Zhang's normalization convention 
for Sabcd] see Eqs. ([2?7|. 



£°a = 'yaScdf^O,'^^'^ 

£1^ = 27„'=7,''f,d,n'= + -yabE" 



£l = 27,'=f,de/n'^n'=n^ 

£'ab = ^-fa^lb^Scdefn-^f + -f^bS' 



TABLE II: Irreducible tidal potentials: type-Z?, odd-parity. 
The superscripts q, o, and h stand for "quadrupole," "oc- 
tupole," and "hexadecapole," respectively. The tidal mo- 
ments Bab, Babe, and Babcd are the STF tensors defined in the 
text. The factors of | and ^ are inserted to respect Zhang's 
normalization convention for Babe and Babcd', see Eqs. (12. 61) 



and (2.71 







Blb^ 


eap,^l^B\-i% + ebp^n^B^.y^a 


131 = 


S^apq^^B'^^^Q, n 


Blb = 


. Ueapgn'^B'^yt + ebpqQ.m\e-i\)Q,^ 


B'a = 




B'ab = 


- fitap^n^B'^^^Y, + ebp^n^B^^fj^aWnf 



TABLE IV: Irreducible tidal potentials: type-Z3Z3, even- 
parity. The superscripts m, q, and h stand for "monopole," 
"quadrupole," and "hexadecapole," respectively. 



BpqBP" 



Q" = Bp^.B^'^^Q'^n" 
Ql=laBp(,B''^^n'' 

Qlb=2la^b''Bp^eB\,+^abQ'' 



Q" = B(edBef)n-'Q,'^Q.'nf 

&a=la'B(,iB,f)n''n^^f 

Qlb = 2^a''lb''B^edBef)n-nf +^abQ' 



TABLE V: Irreducible tidal potentials: type-£S, even-parity. 
The superscripts d, and o stand for "dipole" and "octupole," 
respectively. 







Qi = 




9" = 


epq{.£''iB\^n''n''n'' 


Ql = 


^a'ep,^,£^^B\^^''n^ 


G°ab = 


--'2la^^b%qic£^dB\^^'+'yab Q" 



and adopt TZ as the single scale associated with the tidal 
environment. In this sloppy notation we write 



^ab ^ Bab 



1 

TV' 



^abc ^ ^abc ^ ^ab ^ ^ab 



1 

TV' 



(2.12a) 
(2.12b) 



When the motion is slow we have that C <^TZ, and there 
is a wide separation between the two length scales. When 
the motion is relativistic, however, all scales are compa- 
rable to each other. 

To simplify our notation in later portions of the paper, 
we choose to eliminate the distinction between the dif- 
ferent tidal scales. We therefore set V ^ I, £ ^ T ^ TZ, 



TABLE III: Irreducible tidal potentials: type-ff, even- 
parity. The superscripts m, q, and h stand for "monopole," 
"quadrupole," and "hexadecapole," respectively. 

-pm _ c cpq 



•pq 

-pq 
' a 



- £ p ^ c £ ^ ^ 

' !a ^p{c^ d) 



r:, = 2ja^'y,%^,£^^^+'yabV' 

V =£i^cd£ef)^''^''^''^^ 



-ph - 
' a 

' ab 



-ia'£^,i£,f)^''n^nf 

--2'^a'lb''£{cd£el)^^^^ +labV'' 



Eabca - Babcd - Sabc ^ B^bc Sab ^ B^b ^ (2.12c) 



We emphasize that this relaxation of our notation is sim- 
ply to save writing. For example, an error term that 
should be written 0{r^Tl^^ C^^T^^) will be condensed 
to the simpler expression 0{r^ /TV). The form of the 
equations will always allow us to determine the order of 
magnitude of each term in relation to the complete set 
of scaling quantities. 



TABLE VI: Irreducible tidal potentials: type-i£.B, odd-parity. 
The superscripts q, and h stand for "quadrupole" and "hex- 
adecapole," respectively. 

ni = eJ"'np£,^,B\~^^l^ 

Til, - ea^'>np£,^,B^^)-i% + er^lp£,^,B\~^-i\ 

Ti)i = eJ"^Q.p£ (^qcBde)^'^^'^^'^ 

H\ = {ea^''np£(^dBcf)-i\ + e,^''np£(,,B,f)-i^aWnf 
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C. Parity rules 

In the context of this work, a parity transformation is a 
change of tetrad vectors described by — > u" and e" — > 
— e"; the transformation keeps the permutation symbol 
unchanged: eabc — > £abc- Under the transformation the 
frame components of the Weyl tensor change according 
to CaobQ Caobo and Cabco -Cabco- Wc also have 



C, 



aObO\c 



-c, 



C'aO&Ojcdi C^abcO\de 



a060|ci C'af)cO|d 



c, 



-a 



abcO\d: 



and C, 



abcO\de ■ 



From Eqs. (2.4 1, (2.6 I 



aObO\cd 



and (2.7 1 we deduce the transformation rules 



Sab ' 
Bab 



Sab J Sabc ^ 
' Bab J Babe 



Sabc7 
Babe, 



Sabcd 
Babcd 



Sabcd, 
—Babcd 



(2.13a) 
(2.13b) 



for the tidal moments. We say that Sabi Sabc, and Eabcd 
have even parity, because they transform as ordinary 
Cartesian tensors under a parity transformation. And we 
say that Bab, Babe, and Babcd have odd parity, because 
they transform as pseudotensors. 



D. Tidal potentials: Cartesian coordinates 

For the purposes of writing down the black-hole metric 
in Sec. |III[ we involve the tidal moments in the construc- 
tion of tidal potentials that will form the main building 
blocks for the metric. We first achieve this with the help 
of a system of Cartesian coordinates x° that we assume 
is at our disposal. In the following subsection we convert 
the potentials to spherical coordinates {r,d,<j>). 

We introduce 



:= x°-lr 



(2.14) 



as the radial unit vector, with r := y 6abX°"^ denoting 
the usual Euclidean distance. We refer to the radial direc- 
tion as the longitudinal direction, and to the orthogonal 
space as the transverse directions. 

We wish to combine the tidal moments with il'^ so as 
to form scalar, vector, and (rank-two, symmetric) tensor 
potentials that satisfy the following properties: 

1. Each potential is an element of an irreducible repre- 
sentation of the rotation group labeled by multipole 
order /. 

2. Each scalar potential transforms as such under a 
parity transformation. 

3. Each vector potential transforms as such under a 
parity transformation, and is purely transverse, in 
the sense of being orthogonal to il'^. 

4. Each tensor potential transforms as such under a 
parity transformation, and is transverse-tracefree, 
in the sense of being orthogonal to 51" and having 
a vanishing trace. 



TABLE VII: Spherical-harmonic functions y'™ and harmonic 
components Am involved in the decomposition of .A*^'' := 
Ak^.-ki^''^ ■■■fl''' = T.m^^^Y'^"'- The functions are real, 
and they are listed for the relevant modes I — 1 (dipole), 1 — 2, 
(quadrupole), 1 = 3 (octupole), and I = 4 (hexadecapole) . 
The abstract index m describes the dependence of these func- 
tions on the angle (j>; for example y'--^^ is proportional to 
sin2<^. To simplify the expressions we write C := cos 9 and 
S ~ sin 6*. The harmonic components are expressed in terms 
of the independent components of the STF tensor Aki^-'ki ■ 








= c 


A'' - 


= -43 




Ic 


— S COS 




= Ai 


yi 


Is 


— S sin 




= Al 


y2 





= 1 - 3C^ 


Al -. 


= i(-4ll +-422) 


Y2 


Ic 


= 2SC cos<p 




= -4l3 


Y2 


Is 


= 25Csin<^ 


A. 


= -423 


Y2 


2c 


= cos 20 




= i(-4ll --422) 


Y2 


2s 


= sin 20 


A. 


= -4l2 


Y3 





= C(3 - 5C^) 


AI-. 


= 3(-4ll3 + -4223) 


y3 


Ic 


= |S(1 - 5C^)cos0 


-4ic 


= 5(-4lll +-4l22) 


y3 


Is 


= 15(1 - 5C^)sin0 


-4is 


= 3(-4ll2 + -4222) 


Y3 


2c 


= 35^0 008 20 


-4°2c 


— 2 (-4ll3 — -4223) 


Y3 


2s 


= 3S^Csin20 


-4°2s 


= -4l23 


y3 


3c 


— cos 30 


A. 


= i(-4lll - 3-4l22) 


y3 


3s 


= sin 30 




= j(3-4ll2 - -4222) 


y4 





= i(3 - 30C^ + 35C*) 


A^ - 


- \ (-4llll + 2^1122 + -42222) 


y4 


Ic 


= 25C(3- 7C^)cos0 


A. 


— ^(-4lll3 +-4l223) 


y4 


Is 


= 2SC(3 - 7C^)sin0 


As 


— ^ (-4ll23 + -42223) 


y4 


2c 


= S^(l - 7C^) cos 20 


^2c 


— ^(-4llll — -42222) 


y4 


2s 


= S^(l - 7C^)sin20 


A'' 
-^2s 


— -4lll2 + -4i222 


y4 


3c 


= 45^ C cos 30 


A^ 


— ^(-4lll3 — 3-4l223) 


y4 


3s 


= 4S^Csin30 


As 


— ^(3-4ll23 — -42223) 


y4 


4c 


= cos 40 


a'' 


— 1 (^1111 — 6-4ll22 H- -42222) 


y4 


4s 


= S'^ sin 40 


A^ 

^4s 


— ^(^1112 — -4l222) 



The first property implies that each potential will sat- 
isfy an appropriate eigenvalue equation that depends on 
the multipole order and the tensorial rank of the po- 
tential; these are displayed in Eqs. (A3 1 and (AlO) of 



Appendix |A] Under a parity transformation the longitu- 
dinal vector transforms as fi" — fi", a scalar potential 
remains invariant, a vector potential changes sign, and a 
tensor potential remains invariant. To aid the construc- 
tion of the potentials we introduce 



^\ := 5\ - n'^nb 



(2.15) 



as a projector to the transverse space orthogonal to fi"; 
this transforms as jab lab under a parity transforma- 
tion. 

The required potentials are displayed in a number of 
tables. In Table|l]we list the tidal potentials that are con- 
structed from the even-parity tidal moments £ab, £abc, 
and Eabcd- In Table [II] we have the tidal potentials that 
are constructed from the odd-parity tidal moments Bab, 
Babe, and Babcd- In Table [Hi] we list the potentials that 
arise from the bilinear combination Sabred of even-parity 



moments. In Tabic IV wc have the potentials that arise 



9 



TABLE VIII: The first column lists the harmonic components 
of type-f, even-parity tidal potentials, as defined in Table |l] 
The second column lists their expansions in scalar, vector, 
and tensor harmonics. The components of £^ come with an 
additional factor of 2 to accommodate Zhang's choice of nor- 
malization; see Table |l] 



c^ _ 
"^0 - 


= +£22) 






c-q 


= £13 


f = 


cq \r2,m 




= £23 


cq _ 

- 


1 ^ cq 
" 2 Z^m ^ A 




= 1 11 — £22) 


cq 

'-'AB 


\ ' cq T/'^jm 

^ Z^m '-■m^ AB 




= £12 






£■1 -- 


- 2 (£^113 + £223) 






£-ic 


= 1 (fill + £122) 






£-\s 


= 1 (£112 + £222) 


£° = 


CO \^3,m 
Z^m '-'m-' 




= ^{£113 — £223) 


S'a = 


1 \ c'o 'v/'3,Tn, 
3 Z^m ^m'^ A 


£-2s 


= £123 


£-AB 


1 CO \^3,m 

— 3 Z^m ''m-' AB 


£3c 


= \ (£111 — 3i£l22) 






F?, 
"--38 


~ 4(3^^112 ^^222) 






C-h . 

•^0 - 


= \ (£1111 + 2^1122 + £2222) 






"--Ic 


= £1113 + f 1223 






ch 


= f 1123 + f 2223 






fh 

02c 


= £1111 — £2222 


£' = 


ch "^^4,771 

Z-(m J 


ch 
<^2s 


= 2(f 1112 + f 1222) 


ch _ 

C-A - 


1 ch T^4,m 
4 Z^m A 


oh 

^3c 


= |(iSlll3 — 3^1223) 


ch 

^AB 


1 Y^ ch \^4,m 

~ 6 Z^m ^rn'- AB 


C-3s 


= |(3i£ii23 — £2223) 






oh 
^•4c 


= 1 (full — 6^1122 + f 2222) 






oh 
•--48 


= £1112 — £1222 







TABLE IX: The first column lists the harmonic components 
of type-Z3, odd-parity tidal potentials, as defined in Table [H) 
The second column lists their expansions in scalar, vector, 
and tensor harmonics. The components of B° come with an 
additional factor of | to accommodate Zhang's choice of nor- 
malization, and those of B'' come with an additional factor 
f; see Table |ll] 





= |(Sll +Z322) 








= Bis 


B" = 


\ ^ Kq ■\/'2,m 




= B23 


B\ = 


1 ^ — > nn f .7-2. 771 

- 2 Z^m '-'m^A 


62c 


= |(Sll - Z322) 


Bab 


Y"^ loq v"2,m 
— l^m '-'m^AB 




= B12 






B°o -- 


= 1 (S113 + S223) 






Blc 


= |(Siii +S122) 






Bis 


= |(Sll2 + S222) 


B° = 


\ — ^ via T 7T), 


B°2c 


= 3 (S113 — S223) 


B°A = 


1 \ fJO \^ 3 , ?TT- 

= 3 Z^m "m-^A 


132s 


4 f2 

- 30123 


Bab 


= 3 "m-^AS 


Bic 


= |(Z3iii - 3^122) 






Bis 


= |(3Sii2 — S222) 






B'o-- 


= |(Biiii + 2B1122 + B2222) 






B\c 


= 1 (S1113 + S1223) 






B'ls 


~ |(^3ll23 + S2223) 






B'2c 


= 1 (Bull — Z32222) 


B' = 


Y^ f2h T^4,m 
Z^m 


Bis 


— ^ (S1112 + B1222) 




1 Y^ »3h v^'^ 
- 4 Z^m "m^A 


Bl^ 


= §(^1113 ^ 3Z3i223) 


B'ab 


1 Y^ f2h Y"4.7n 
^ 6 Z^m "m^AS 


Bis 


— §(3Z3ii23 — ^32223) 






Bt 


= Y^(Siiii — 6S1122 + S2222) 






Bis 


= §(^1112 — B1222) 







from the bilinear combination Babied of odd-parity mo- 
ments. In Table [V| we list potentials that arise when we 
combine Sab and Bab into an even-parity structure. And 
finally, in Table |VI| we combine them into an odd-parity 
structure. 

We illustrate the construction of the potentials by ex- 
amining a few examples. The simplest cases involve the 
tidal quadrupole moment Sabi which transforms as a ten- 
sor under a parity transformation. The associated scalar 
potential is £ab^°'^^', this transforms appropriately un- 
der a parity transformation, and because £ab is tracefree 
this satisfies the scalar eigenvalue equation with 1 = 2. 
To get a vector potential we first form Ecb^'', which has 
the appropriate number of free indices. We next multiply 
this by j'^a to make the vector transverse, and verify that 
the final object satisfies the vectorial eigenvalue equa- 
tion with 1 = 2. To get a transverse tensor potential we 
first form 'y'^al'^b^cd, and we remove its trace by adding 
^Jabi^cd^"^''')', the result satisfies the appropriate eigen- 
value equation with I = 2. The method generalizes easily 
to octupole and hexadecapole potentials, and the final 
results are displayed in Table |Tj 

We next examine the potentials associated with the 
odd-parity quadrupole moment Bab- The combination 
Bab^"^^^ transforms as a pseudoscalar under a parity 
transformation and does not, therefore, satisfy the crite- 



TABLE X: Harmonic components of type-££, even-parity 
tidal potentials, as defined in Table |III| The spherical- 
harmonic decompositions are V = X^mP™^^ '"' ^a ~ 

1 Y^ 'pq y2,m ^q Y^ y-2,m ^h Y^ -ph y4,m 

2 Z^m ' rn-' A ' ' AS Z^m ' m-' AB ' ' Z^m ' ^ ' 

n = \ j:,^rLY^'"', and p'Xb = | j:,^riY^'B ■ 



V- = &{£^,f + 2{£l^f + 2(g;j^ + 2{£l^f + 2{£l^f 



-pq . 


= -(£^5)2 - ii£ir - msf + wicf + u^is 


r 


-pq 

' Ic 


cq cq _i_ cq cq _i_ cq cq 




-pq 

' Is 


cq cq _i_ cq cq cq cq 






= 2£^o£L + Wlcf~'2(£lsf 




vis 






Po = 


= m^'of - mi)' M^is? + M^icf + is 


{£ls? 


Pjc 


6 cq cq _i_ 1 cq cq i 1 cq cq 

~ 7*-'0*-'lc "1" 7*-'lc*-'2c "T" 7*-'ls*-'2s 




vis 


6 cq cq _i_ 1 cq cq i cq cq 

~ 7*-'0*-'ls ~r 7*-'lc*-'2s 7*-'ls*-'2c 




VI, 


= ff^lf^, -!(£:? J' + f(£:;j2 




vis 


6 cq cq 4 cq cq 

~ 7*-'0*-'2s 7*-'lc*^ls 




VI 


1 cq cq 1 cq cq 

~ 2*-'lc*-'2c 2*^ls<-'2s 




VI 


1 cq cq 1 1 cq cq 

~ 2*-'lc*-'2s ^ 2*^ls*^2c 




P4c 


= '2(£lc? - U^lsf 




VI 


cq cq 

— *--2c*--2s 
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TABLE XI: Harmonic components of type-BB, even-parity 
tidal potentials, as defined in Table |IV| The spherical- 
harmonic decompositions are Q*" = QlnY^'"", Q'J^ = 



and Q^B = |E™Q5.y4T 



or = 6(SS)2 + 2(S1J2 + 1{B\,f + 2(SL)' + 2(SL 



= 2Z?seL + e^SL 



-'2s 



ACRl ^2 _ ^ 



4K7q m 
t'-'Ic'-'Is 



= \B\j3i, + i^LB: 



2c 
2 



ria that were formulated previously. To obtain a suitable 
potential we must involve the permutation symbol, and 
the simplest allowed combination is Capq^^B'^^il'^. This 
transforms as a vector under a parity transformation and 
satisfies the appropriate eigenvalue equation with 1 = 2; 
we therefore include it as one of our building blocks. To 
form a tensor potential we remove fl'^, replace it with 
the projector and symmetrize the indices; the end 
result satisfies all the properties required of a tensor po- 
tential. The method generalizes easily to octupole and 
hexadecapole potentials, and these are listed in Table [iT] 
The bilinear potentials are constructed by forming STF 
products of the quadrupole moments Sab and Bab and 
combining these with appropriate factors of fi", 7";,, and 
tabc- For example, the STF products that can be formed 
from two factors of £ab are the scalar £pq£P'^, the rank- 
two tensor Spi^a^^b) ' ^^'^ rank-four tensor £{ab^cd) the 
associated potentials are listed in Table and those 
of Table |IV| follow from similar manipulations. When 
Sab and Bab are both involved we must be mindful of 
the parity rules; for example, the hexadecapole poten- 



tial E/^abBcd)^"'^''^'^^''^ is ruled out because it transforms 
as a pseudoscalar under a parity transformation. The 
appropriate combinations must involve the permutation 
symbol, and these are displayed in Tables [V] and |VI| 



E. Tidal potentials: Angular coordinates 

A transformation from Cartesian coordinates to 
spherical coordinates (r, 9, (/>) is effected by 



ai-nA\ 



(2.16) 



in which the longitudinal vector is now parameterized 
by two polar angles 9^ = {9,(j)). Explicitly, we have that 
il"' = [sin cos 0, sin sin (/), cos 0] . The transformation 
imphes that dx°- /dr = fi" and dx°- /d9^ 



rfi^, with 



We note the useful identities 



(2.17) 



(2.18a) 
(2.18b) 
(2.18c) 



Here ^ab = diag[l, sin^ 9] is the metric on the unit two- 
sphere, and is its inverse. We introduce Da as 
the covariant-derivative operator compatible with ^Iab, 
and eAB as the Levi-Civita tensor on the unit two-sphere 
(with nonvanishing components ee^ = —ecf,0 = sin0). We 
adopt the convention that upper-case Latin indices are 
raised and lowered with and ^ab, respectively. Fi- 
nally, we note that Dc^ab = Dqcab = 0. 

We convert the vector and tensor potentials from their 
initial Cartesian incarnations to angular-coordinate ver- 
sions by making use of the transformation matrix f2^. 
We thus define 



CI — PiO" 
^A ■— ''a^Ai 



CI — Ct Qa Qb 



(2.19) 



and apply the same rule to all other potentials; for ex- 
ample, G°AB '■— ^ab^A^s- After this conversion the tidal 
potentials become scalar, vector, and tensor fields on 
the unit two-sphere, and they depend on the angular 
coordinates 9^ only. It is easy to show that the con- 
version can be undone; for example E"^ = f^fi^, with 

^a '■— ^ab^^^^'s- 

The tidal potentials can all be expressed in terms of 
(scalar, vector, and tensor) spherical harmonics. Let y'™ 
be real-valued spherical-harmonic functions (as defined in 
Table |VII[ ) . The relevant vectorial and tensorial harmon- 
ics of even parity are 



yj" DaY^" 



AB • — 



1 



DaDb + -l{l + l)nAB 



(2.20a) 
(2.20b) 



notice that ^^^Y^JJ^ = by virtue of the eigenvalue equa- 
tion satisfied by the spherical harmonics. The vectorial 
and tensorial harmonics of odd parity are 



-^A 



-e/DgY 



Im 



(2.21a) 



^AB ■~ ~2 



Da]DcY'"' = 0; (2.21b) 



the tensorial harmonics also are tracefree: 

flABj^im^ = 0. The decomposition of the tidal poten- 
tials in spherical harmonics is presented in Tables |VIII| 
[TX) [Xj [XT] |XIT] and panl A derivation of these resuhs is 
presented in Appendix |A[ 
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TABLE XII: Harmonic components of type-SB, even-parity tidal potentials, as defined in Table |V] The spherical-harmonic 
decompositions are G" = GlnY^'"', Q\ = E,„ C^i'", Q° = E™ Smi"'''", Qa = | E„ e^i^'", and ^ab = | E„ Sm^lr- 



= (3^3 - f^jsL - - f;,(35s - + s^bi^ 
= -(3gg + gpglc + g;c(3gS + + g^^L - g^^L 

/^o 2 cq K?q _i_ 2 cq K?q _i_ 1 cq Rq l cq K?q 

= -Igf-JSL - IgJc^L + \£lA-iBi + Bi:) - l£l,B\, + \£IM. 
Gls = f g^e^ - |g;,(2SS - BIJ + igJ.SL - |g^,Z?L - |gL6L 

= £o^2s — g^^O 

gL = -g,ie^. + gLsS 

= -|g;.SL - Ig^s^, + ig^eL + Ig^^L 

/^o Ic^iqq Icq f2q 1 cq K?q _i_ 1 cq K?q 

t/3s — 2C-lc'-'2c 2'-'ls'-'2s 2'-'2c'-'lc "r 2'^2s'-'ls 



TABLE XIII: Harmonic components of type-gfi, odd-parity tidal potentials, as defined in Table |VI[ The spherical-harmonic 
decompositions are = | Em'^m^l''", Ti-^g = Em^™^As . = 1 Em^m^l'"', and Hab = | Em^m^l's • 



oyq >^q^?q Icq K?q Icq K?q i 1 cq K?q i 1 cq K?q 

''■0 ^ '-'O'^O 6 Ic'-^lc 6 Is'-^ls I 3'--2c'^2c ' 3'-'2s"23 

T~(-lc = ^ 2^0^1c ^ 2^1c(^0 ^ ^2c) + 2^1s^2s + 2^2c^lc + 2^2s^ls 

= -|g,^,SL + ig;,eL - |gL(sS + ^L) - Ig^cSL + ^L^lc 
HI, = g^^3L + |g?cSL - IgJ.fiL + g^eS 

'^2s = gp^L + 2^1c^ls + 2^1a^le + gL^O 

■j-Zh IScqj^q 4 cq K?q 4 cq »2q i 1 cq f^q i 1 cq f^q 

''•0 ^ 35''0'-'0 ~ 35'^lc^lc ~ as^'ls'-'la + 35'^2c'-'2c + 35''2s'-'2s 

T~(.\c = f go^lc + nglc(6So + ^2c) + ngls^L + ng2c^lc + I4g2s^ls 

= fg,^,eL + h^l.^ls + h^Lm - - rjg^cCL + ^jg^^L 

'i_/h 3 oq Kjq 2 r2^ i ^ cq K?q i 3 cq K?q 

n.2c — 7^0^2c ~ T^lc^lc ^ T^ls^ls ^ 7^2c^0 

o-zh 3 cq K?q 2 cq K?q 2 cq K?q _i_ 3 cq K?q 

n-2s — 7^-0^23 ~ 7^1c^ls ^ 7^1s^lc ^" 7^2s^0 

^3c = ig^^L - ig^CL + ig^e^L - igL^L 

^3. = igJ.SL + ISIBI, + \£l,Bl^ + \£l,B\, 

-Tjh 1 cq K?q 1 cq K?q 

'l4c — 2'^2c*-'2c 2'-'2s'-'2s 

'Tjh 1 cq K?q I Icq Kq 

'1-4S — 2'-'2c'^2s "T 2'-'2s'^2c 



III. GEOMETRY OF A DEFORMED BLACK 
HOLE 

We construct the metric of a tidally deformed black 
hole in two steps. In the first step we continue to think 
of a smooth timelike geodesic 7 in a vacuum region of an 
arbitrary spacetime, and we construct the metric of this 
(background) spacetime in a neighborhood of the world 
line. We denote this neighborhood by A/", and we require 
that it be small in comparison with the length scale TZ 
that characterizes the tidal environment: we demand that 



r < 7^, 



(3.1) 



where r is a radial coordinate (to be introduced below) 
that measures distance to the world line. In the second 
step we insert a black hole of mass M into the background 
spacetime, and place it on the world line 7. For the con- 
struction to be successful it is necessary that the world 
tube traced by the black hole fit comfortably within A/", 
and this can be achieved when 



This condition implies both that the black hole is weakly 
perturbed by the tidal environment, and that the world 
tube is small when viewed on the scale TZ of the external 
spacetime. In this situation it makes (approximate) sense 
to say that the black hole moves on a world line 7. 
The metric of the background spacetime is presented 



in Sec. Ill B and the metric of the black-hole spacetime is 



(3.2) 



presented in Sec. |III C| These subsections contain what is 
truly a presentation of the metrics, and there the reader 
will find no trace of a derivation of our results. The cal- 
culations that lead to the metrics are quite lengthy, and 
derivations are relegated to Secs.|V)|Vl) and Appendix [D| 
We begin in Sec. |III A| with a description of our coordi- 
nate systems. 



A. Light-cone coordinates 

We work with a system (w, r, 6, (f>) that is specifically 
tailored to describe the geometry of light-cone surfaces. 
We refer to these as light-cone coordinates. In the case of 
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the background spacetime (Fig.[lJ, we consider past light 
cones that converge toward the world line 7, so that the 
apex of each cone coincides with a point on the world line. 
In M the light cones provide spacetime with a foliation 
by null hypersurfaces, and each light cone is generated by 
a congruence of null geodesies. (We assume that caustics 
do not develop within M , except at the apex of each light 
cone.) The coordinates are intimately tied to the light 
cones and their generators: 

1. The advanced-time coordinate v is constant on each 
light cone, and its value on a given light cone is 
equal to proper time t on the corresponding point 
of the world line. 

2. The angular coordinates 9^ — {9, (f) are constant 
on the null generators of each light cone; the angles 
refer to a set of axes that are aligned with the basis 



vectors e" introduced in Sec. Ill A 



3. The radial coordinate r is an afhne parameter on 
the null generators of each light cone, normalized in 
such a way that the metric takes a Minkowski form 
in the immediate vicinity of the world line; note 
that r decreases to zero as the generators converge 
toward the world line. 

To illustrate the meaning of the coordinates we examine 
the simple case of an observer at rest in flat spacetime. 
The observer's proper time r is equal to coordinate time 

the advanced-time coordinate is obtained by the trans- 
formation V = t+r, and (r, 9, <j)) are obtained in the usual 
way from the Lorcntzian coordinates x". The metric is 
ds^ = -dv'^ + 2dvdr + r'^dn'^, with d^"^ := d9'^ + sin^ 9 dcj? 
denoting the metric on the unit two-sphere. It is easy 
to verify that the light-cone coordinates satisfies all the 
properties listed previously. 

In the case of the black hole spacetime (Fig. [2]), the 
light cones no longer converge toward a world line. In- 
stead they converge toward the horizon, which traces a 
world tube in spacetime. They still, however, provide M 
with a foliation by null hypersurfaces, and each light cone 
is still generated by a congruence of null geodesies. The 
light-cone coordinates keep most of the properties listed 
previously: 

1. The advanced-time coordinate v is constant on each 
light cone. 

2. The angular coordinates 9^ = {9, ip) are constant 
on the null generators of each light cone. 

3. The radial coordinate r is an affine parameter on 
the null generators of each light cone. 

Since the calibration by the world line is no longer avail- 
able, the coordinates have lost some of their rigidity. We 
can, however, restore most of this rigidity by imagin- 
ing two spacetime foliated by past light cones. The first 
is the background spacetime, with its fully specified set 
of light-cone coordinates. The second is the black-hole 



V = const 




FIG. 1: Light-cone coordinates centered on a world line 7. 
The figure shows a light cone v = constant that intersects 
the world line at the point z{v). It shows also one of the 
light cone's generator, along which 6^ is constant; the afline 
parameter r decreases to zero as the generator approaches the 
world line. 



V = const 




FIG. 2: Light-cone coordinates centered on a black hole. The 
figure shows the world tube traced by the black-hole horizon. 
It shows also the light cone v = constant and one of its gener- 
ators. As in Fig. 1, the angles 6^ are constant and the afline 
parameter r decreases on the generator. 



spacetime, with its own set of light-cone coordinates. 
Far from the black hole, where r is much larger than 
M but still much smaller than TZ, the gravitational in- 
fluence of the black hole is small, and light rays behave 
there just as they do in the background spacetime. We 
can therefore tune the black-hole coordinates so that the 
asymptotic description of the null generators agree with 
the background description. This correspondence does 
not permit the complete specification of the coordinates; 
there remains a limited freedom to redefine the coordi- 
nates without changing their geometrical meaning. We 
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have thoroughly exploited this freedom to simplify the 
form of the black-hole metric. 

The light-cone coordinates are best described in terms 
of the radial distance r and the angles 9^. It is use- 
ful, however, to introduce also a variant of the light-cone 
coordinates that involves the quasi-Cartesian system a:" 
instead of the quasi-spherical system {r,6^). The coor- 
dinates a;° are defined in the usual way by Eq. (2.16 1. In 



the flat-spacetime example examined previously, we find 
that the components of the Minkowski metric are given 
by 77«D = -1- Vva = ^a, and rjab = Jab when presented in 
the light-cone coordinates (w, a;"). It should be noted that 
the metric is mildly singular at x'^ — 0, because its value 
depends on the ambiguous direction of the longitudinal 



vector fi". This singularity persists in the background 
spacetime, and dealing with it requires some care; but it 
is not a serious obstacle to any computation. 



B. Background spacetime 

The metric of a vacuum region of spacetime that sur- 
rounds a timelike geodesic 7 can be presented in the 
light-cone coordinates introduced in the preced- 

ing subsection and expressed as an expansion in powers 
of r/TZ ^1. It is given by 



3 



21 



2 

15 



An'' 



9va — 



1 



-r 



7 



3. r 



1 

6 
Q 



15 



12 

-0(5), 



(3.3a) 



75 

gab = lab ^ 
1 



K + ^r^ni + ^r^Vl + llQl) + \t^GI - ^^'(^ + Sa) + 0(5), (3.3b) 



21 

\r^[^lb-Blb) + Y/[^lb 
^,r^Kb~B^ab) + i-,rSab{V 
'^r'n:,-l^r^labG^+lr^ 



Bib) - I 



10 1 R 



B 



'abl 



20 



y'lab{V'+Q')+0{b) 



^^r'lab G' - ^^r^lab{V^ + Q') - ^^r'iKb - Ql 



■abj 



(3.3c) 



The metric features the tidal potentials encountered in 
Sec. II D and the notation 0(5) indicates that the error 
terms are of order (r/TZ)^. The world line is situated at 
x° = 0, and V is proper time on 7. 

The terms of order (r/TZ)'^ involve the quadrupolc tidal 
moments Sab and Bab, which are now expressed as func- 
tions of advanced time v instead of proper time r. At 
order (r/TZ)^ we find terms involving Sab and Bab, the 
time derivatives of the quadrupole moments, as well as 
the octupolc tidal moments Sabc and Babe- And at or- 
der {r/TZ)^ we see occurrences of Sab, Bab, Sabc, Babe, 
the hexadecapole tidal moments Sabcd, Babcd, and bilin- 
ear combinations of the quadrupole moments. We see 
that the metric involves the majority of the tidal poten- 
tials listed in Tables |l] |lT] |llT] [TV] [V] and [VT] Exceptions 



are V\ H^, Gab, Kb, Qab^ and W,,; these could have 
occurred in the metric, but they happen to be ruled out 



by the Einstein field equations. 



A useful way to view the metric of Eqs. (3.3) is to 



recognize that it possesses the correct number of freely- 
specifiable functions of time to describe correctly, and in 
sufficient generality, the geometry of a vacuum region of 
spacetime in a neighborhood A/" of a timelike geodesic. 
The functions of time are contained in the tidal mo- 
ments; we have 10 of them in Sab and Bab, 14 in Sabc 



and Bnhc, and 18 in 



bed 



and Br 



bed' 



The total is 42, the 



correct number for a metric constructed through order 
(r/TZ)'^. These functions encode meaningful information 
about the spacetime geometry; as we saw in Sec. |II A| 
they represent components of the Weyl tensor and its 
derivatives evaluated on the world line 7. 

In quasi-spherical coordinates {v,r,6^) the nonvanish- 
ing components of the metric are 



9vv 



-1 - r^S 
2 

' 15 



-r^S" 



3 21 
1 



-r^S° - — 
6 12 



4 oh 



15 ^ 



r^G' + -r^Q" + -r^G° - -r''{V' + Q') + 0(5) 



(3.4a) 
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gAB = r Hab 
1 
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21 
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AB ^AB 



18 

ABj-r^''^AB{V-+Q 

2 



r5(n + QA)+^0(5), 



1 

15 



B\ 



15 

AB ~ ^°AB 



32 16 



AB ^Ab) 



3 

20 



rH^AB 



(3.4b) 
(3.4c) 



„6 



45 



nAB{V' + Q' 



^0(5). 



(3.4d) 



This can be obtained from Eq. (3.3 1 by applying the 
transformation rules described in Sec. Ill El The state- 
ments that = 1 and QtA — are exact, and they 
follow from the light-cone nature of the coordinate sys- 
tem. Here the metric features the tidal potentials listed 
in Tables [VlIIl |IX] [Xj |Xll |XlT] and pOlIl When deahng 
with the angular coordinates it is convenient to rely on 
the decomposition of the potentials in scalar, vector, and 
tensor harmonics. 



C. Black- hole spacetime 



The metric of a tidally deformed black hole resembles 
closely the background metric presented in the preceding 
subsection. In the quasi-Cartesian coordinates (u, x"") it 
is given by (/ := 1 — 2M/r) 



+ '^r^9{Q' + '^r^plV + qlQ^) + ^rVi^" - \r^{p\V' + q\&) + 0(5), (3.5a) 
9va ^^a~ \r\e\El ~ h\Bl) + \r\elEl - blBl) - Jr3(e^f; - - ^r^elE^ - blBl) 
+ Ir^elE:, blBl) ^r^e^E^, - b\B^,) - ^,4^.^. + ^^r^Km 

+ ^^^\Pl'Pl + + \r^92Q°a + r^h\K " ^^'(^2^ + q\Q'a) + 0(5), (3.5b) 
9ab = lab - \r^e'\El, - b}Bl,) + ^r'{elE:, ~ blBl,) ~ ir3(e°£», ~ - \r\elEl, - blB^,) 

01 o 00 

+ ^rSab{v'{P' + + r\v\Vl^ + g^Qj..) + r^h'^\^, + 0(5). (3.5c) 



A comparison with Eqs. (3.3) reveals that there are two 
main differences between the metrics. The first is that 
the black-hole metric involves the radial functions , e° , 



ejj, 6^, 6^, and b\ that are listed in Table XIV as well 



as the radial functions p^, K> C> 9«, .9n> 
and /i° that are listed in Table |XV| The second is that 



the tidal potentials V\U\, Glt^b: 2afc> and H'^,, now 
make an appearance in the black-hole metric. (Recall 
that they were absent in the background metric.) The 
black-hole metric contains the same number of freely- 



specifiable functions of time as the background metric; 
these are contained in the tidal moments Eat, Bab, Eabc, 
Babe, Eabcd, and Babcd, which all depend on advanced time 

V. 

The metric of Eqs. (3.51 contains two fundamental 



length scales, the black-hole mass M and the tidal ra- 
dius TZ. It is assumed that M <C 7?., and the metric is 
valid when r ^TZ. When M ^ the metric becomes the 
background metric of Eqs. ( 
the fact that most of the rac 



3.3 1. This can be seen from 



ial functions approach unity 
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TABLE XIV: Radial functions: linear potentials. The radial functions are expressed in terms of x := r/{2M) and / := 1 — 1/x. 

The dilogarithm function is defined as diloga:: = dt log(t)/(l — t), with logo; denoting the natural logarithm. All radial 
functions vanish at r = 2M, except for = |, = jq, = fey = — |, b°i = —j^, 67 = — 

e\ = f 

4 = / [I + 3^(5 + 12 log x) - j^(27+ 121og.T) + + 

4 = 1 + 2^; (89 + 84 log x) + (431 + 996 log x - 1680 dilog x) - ^ (315 + 282 log x - 210 dilog x) 

+ 125^^(4183 + 2322 log a; - 1260 dilog - ^ - 
el = f 

4 = /[I + 6^(13+ 121ogx) - ^ - ^ - ^] 

e^6 = 1 + (1 17 + 84 log x) + ^ (929 + 1836 log x - 1680 dilog x) - (597 + 387 log a; - 210 dilog x) 

I 387 I 809 I 809 
80x* 480x^ 960a:« 

e^s = 1 + ^(4 + Slog a:) - - (7 + 3 log a:) + ^ 

4 = 1 + 2ik (739 + 420 log x) + ^ (262 + 387 log a; - 210 dilog a;) - (317 + 70 log x) 

- 545^^(1511 + 1161 log X - 630 dilog x) + ^ 

^2 = / [1 + ScbV3 + 60 log x)-^ (479 + 180 log x) + ^ (87 + 20 log x) - ^ + ^] 
e4 = /(l^^) 

65 = / [I + 3^(103 + 60 log x) - ^ (233 + 60 log x) + ^ - ^ - ^] 

e°7 = / + inb^ 

e°8 = 1 + 53^ (85 + 601ogx) - ^(10 + Slogx) + + ^(4 + logx) + ^ 
64 = /(I - ^ + lib) 

67 — ^ 3a: ^ 7a;2 42a:4 

bl = f 

bl = /[I + 6^(7+ 121ogx) - ^ - ^ - ^] 

6^ = 1 + (75 + 84 log x) + (649 + 996 log x - 1680 dilog x) - (879 + 564 log x - 420 dilog x) 

I 141 I 223 I 223 
~'~ 40a;'* ~'~ 160x5 320x6 

^7 = 1-^ 

''8 = 1 + 5^(5 + 61ogx) - ^ - ^(2 + 31ogx) + 

^9 = 1 + 2ik (529 + 420 log x) + 3gi^ (593 + 1128 log x - 840 dilog x) - 3^ (106 + 35 log x) 

- 108^(2357+ 1692 log X - 1260 dilog x) + Jf^ 

= / [1 + ^4S^ (97 + 60 log X) - ^ (227 + 60 log x) + ^ + ^ + ^] 
bl = l+ 5^^(79 + 60 log x) - ^(97 + 301ogx) + ^ + ^(13 + 31ogx) - ^ 

^'4 = /(l - 1^ + 

»,h _ 1 _ A _1_ _5 1_ 

"7 — ^ 3x 7x2 14J.4 



when 2M/r — » 0. Exceptions are p\, g\, p\, q\, h\, and 
h\, which all approach zero; these functions come with 
the tidal potentials that were absent in the background 
metric. In this limit the tidal moments acquire their pre- 
cise relation with the Weyl tensor (and its derivatives) 
evaluated on the world line a;" = 0. 



When M ^ the world line disappears and is replaced 
by the world tube traced by the black-hole horizon. But 

the black-hole metric continues to approach the back- 
ground metric when r is very large compared with M 
(but still much smaller than TZ). In this case the inter- 
pretation of the tidal moments is more subtle. They still 
provide a complete characterization of the tidal environ- 
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TABLE XV: Radial functions: bilinear potentials. The radial functions are expressed in terms of x :— r /{2M) and / := 1 — 1/x. 
All radial functions vanish at r = 2M, except for Pz = = p\ ~ q\ = p\ = ~ j^^, = = — jq, p\ = -g^, cli ~ ^ ifij, and 
/13 = All functions approach unity as a; — > 00, except for p\, p\, q\, gX, h\, and h\, which all approach zero. 
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mcnt, but they are no longer associated with the Weyl 
tensor evaluated on a world line. 

To elucidate the meaning of the tidal moments it is 
helpful to return to the two spacetimes of Sec. |III A} The 
first is the background spacetime with its central world 
line, and the second is the black-hole spacetime in which 
no such world line exists. We assume that the space- 
times share the same set of tidal moments: the functions 
of time contained in 8ab, Bab, £abc, Babe, £abcd, and Babcd 
are the same in each spacetime. An observer in the first 
spacetime can measure the tidal moments anywhere and 
relate them to the Weyl tensor (and its derivatives) eval- 
uated on the central world line. An observer in the sec- 
ond spacetime can still measure the tidal moments, but 
is prevented by the black hole from relating the results 
of his measurements to a world-line Weyl tensor. But 
the second observer is free to imagine that the black hole 
could be removed from the spacetime without altering 
the conditions in the asymptotic region r ^ M, so that 

1 This is a thought experiment that could not be realized physi- 
cally. It is meant to reflect a mathematical procedure in which 
M is taken to zero while keeping the tidal moments fixed. The 
procedure disregards the fact that in most physical applications, 



his measurements could, after all, be related to a world- 
line Weyl tensor.^ This thought experiment provides him 
with a loose interpretation of the tidal moments: They 
can still be related to a world-line Weyl tensor, but the 
Weyl tensor and the world line are fictitious extrapola- 
tions to r = obtained from information available in the 
r ^ M portion of the black-hole spacetime. 

When the tidal perturbation is turned off (by putting 
all the tidal moments to zero, thereby sending TZ off to 
infinity), the black- hole metric becomes 

gvv=~f, 9va = ^a, 9ab = lab, (3.6) 

where, we recall, f — \ — 2M/r. This is the well-known 
Schwarzschild solution expressed in the light-cone coor- 
dinates {v,x°'). 

In quasi-spherical coordinates (w, r, 0^) the nonvanish- 
ing components of the black-hole metric are 

the tidal moments carry a dependence on M that is revealed 
by matching the local black-hole metric to a global metric that 
contains the black hole and the external bodies. 
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+ gSQ^) + r'ip^VAs + qlQAB) + r'h\n\B + r^O{b). 



(3.7d) 



In these coordinates the TZ ^ 00 limit of the metric is 
ds^ — / dv"^ + 2 dvdr + dfi^ , the Eddington-Finkelstein 
form of the Schwarzschild solution. 



IV. GEOMETRY AND DYNAMICS OF THE 
DEFORMED HORIZON 

In this section we extract the consequences of the 
black-hole metric of Eqs. (3.7 1 on the structure and dy- 
namics of the tidally deformed horizon. We begin in 
Sec. |IV A| with a proof that in the coordinate system 
adopted here, the deformed horizon continues to be de- 
scribed by r = 2M. In Sec. |IVB| we display the compo- 
nents of the induced metric on the horizon. In Sec. II V Cl 
we examine the congruence of null geodesies that gener- 
ates the horizon, and derive expressions for its expansion 

we in- 



scalar and shear tensor. And finally, in Sec. IV D 



tegrate Raychaudhuri's equation and calculate the rate 
at which the horizon grows as a result of the tidal inter- 
action. 



Position of the horizon 



The black-hole metric of Eqs. ( 3.7 1 is presented in light- 
cone coordinates {v,r,9^) whose geometrical meaning 
was described in Sec. IIII A| As we saw, the coordinates 



are tied to the behavior of incoming light rays that mesh 
together to form light cones that converge toward the 
black hole. And as we saw, the light-cone coordinates 
are not fully specified; the limited coordinate freedom 
that remains was exploited to simplify the form of the 
metric. 

Concretely, the coordinate freedom was utilized to im- 
pose a set of horizon-locking conditions that are designed 
to keep the black-hole horizon in its usual place. As a 



result, the horizon of a tidally deformed black hole is sit- 
uated at 



= 2M[H-0(5)] 



(4.1) 



in the light-cone coordinates that give rise to the metric 
of Eqs. \3. 71). In Eq. (4.1 1, and in all equations below. 



the symbol 0(5) means that the correction terms are of 
order {M/TVf' . The horizon-locking conditions are 



?«t, = = 5. 



vA 



at r = 2A/, 



(4.2) 



and it is easy to verify that the metric of Eqs. (3.7 1 sat- 
isfies this property. 

The horizon-locking conditions imply that the sur- 
face r = 2M is a null hypersurface in the black-hole 
spacetime. This is seen from the fact that g"^"^ = at 
r = 2 M, w hich follows from the statement that thanks 
to Eq. (4.2 1, — 1 and gAB are the only nonvanishing 



components of the metric at r — 2M. The null hypersur- 
face is generated by a congruence of null geodesies, and in 
Sec. |IV C| we show that the expansion of this congruence 
vanishes. [The expansion scalar is denoted Q, and we 
show more precisely that O = 0{M^ /TZ^). This implies 
that the expansion vanishes through order TZ~'^, which is 
the order of accuracy maintained in the metric and the 
description of the surface.] The surface r = 2M is there- 
fore a stationary null surface, and it is in this sense that 
it describes a black-hole horizon. It is an isolated horizon 
in the sense of Ashtekar and Krishnan [351 IMl ISZ] ■ 

In general we cannot state that r = 2M is an event 
horizon, because the exact position of the event horizon 
depends on the entire future history of the spacetime, 
and is located by tracing light rays backwards in time 
from future null infinity. Because our metric may not 
be accurate for all times, and because it applies only 
to a neighborhood of the black hole, it does not supply 
us with the tools to locate the event horizon. Under 
restrictive assumptions, however, we can go beyond these 
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limitations and prove that r = 2M describes the position 
of the event horizon. 

We consider a situation in which the tidal perturbation 
is switched off after a time v — vi, so that the metric 
returns to its usual Schwarzschild form when v > vi. 
This final Schwarzschild metric possesses the same mass 
parameter M as the original black-hole metric, because 
(as we shall see in Sec. IV D I the total change in mass 
that can result from tidal heating over a time f i '--^ 7?. is of 
order M^/TZ^; this efi'ect is too small to be described by 
our metric. It follows that the event horizon is described 
by r = 2M when v > vi. At earlier times the horizon 
is described by a null hypersurface that joins smoothly 
with r = 2M a,t v — vi. Since r — 2M is null at all times, 



we conclude that the event horizon is always situated at 
r = 2M. 



B. Horizon metric and surface gravity 

The horizon metric is the specialization of the black- 
hole metric to the null hypersurface r = 2M[1 + 0(5)]. 
Because its component along the w-direction vanishes, 
the horizon metric is degenerate and explicitly two- 
dimensional. We denote its nonvanishing components by 
^AB, and these are obtained by inserting r = 2M within 
Eqs. (|3J|. We get 



B' 



32 

45 



AB 



AB} ^ Y^^'^'^i^AB 



Bab 



1 28 

Q-) - —M'^nAB 
' 45 



64 a 
21 



105 ^ 

ab{V + Q') 



AB ^ 

32 



B' 



16 



-M'^nABiV + Q") + —M^V 
1 \ ^ 21 ' 



AB 



160 
"63" 



AB 



208 
"63" 



AB) 



r 



(4.3) 



In Appendix [B] we calculate the determinant 7 of the 
horizon metric. We find the simple result 



7^ = 4M2 sin 6* [1 + 0(5)]; 



(4.4) 



the tidal potentials do not give rise to corrections to the 
horizon's surface element, ^dOdcj). 

The surface gravity k of the perturbed horizon is de- 
fined by the relation k^\7 /^k" = Kk", where the vector 
k" :— [1,0,0,0] is tangent to the horizon's null gener- 
ators. A calculation based on the metric of Eqs. (3.7) 
reveals that 



1 

4M 



1 



16 



32 



M^£,i,c^-n'>n' - —M^SabS'"' + BabB''' 



0(5) 



(4.5) 



We see that the surface gravity is no longer uniform on 
the horizon; the tidal perturbation introduces a variation 
of order (M/TZ)^ over its surface. Notice that this varia- 
tion is associated with changes in the tidal environment. 
When the perturbation is stationary the correction to the 



surface gravity comes from the last term, which is of order 
{M /TZ)* and uniform over the horizon. These observa- 
tions are compatible with the zeroth law of black-hole 
mechanics. 

C. Expansion scalar and shear tensor 

The null generators of the horizon form a congruence 
whose behavior is described by an expansion scalar G 
and a shear tensor aAB', these are defined by (see, for 
example. Sec. Ill of Ref. [3^) 



dvlAB = QlAB + '2.<JAB, 



(4.6) 



together with the requirement that the shear tensor be 
tracefree: 7^^cr^B = 0. The expansion scalar is then 
equal to the trace of dy^AB, so that 6 = \^~^dy^, where 



7 is the metric determinant. Equation (4.4 1 implies that 



e = o{Afin^). 



(4.7) 



This means that up to this level of accuracy, the surface 
r = 2M is foliated by apparent horizons. 



With this we find that Eq. (4.6 1 reduces to gab 
\dv-iAB + 0{MyTV^). This is 



+ —M^nABiV-" + or) - T^M^^riAB Q' - —M'^flABiV' + Q") + -^m^^ab Q" 
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+ hl'^nAB {V" + + I^M^Kb - Q\b + ^^'^'^As + MO{(S) . (4.8) 

I 



D. Tidal heating 

More information about the expansion scalar can be 
obtained by integrating Raychaudhuri's equation (see, for 
example, Sec. Ill of Ref. [39, ) 



AB 



(4.9) 



in which we can neglect — 0{M^^ /Ti}^). The squared 
shear contains terms that begin at order and the 

neglected terms are of order M'' j'R? and higher. The so- 
lution to Eq. (4.9 1 will be an expression for that begins 
at order /TV , and neglects terms of order jV? and 
higher. Given this degree of available accuracy, it is ap- 
propriate to set 



1 

4M 



(4.10) 



in Eq. (4.9 1; the fractional corrections of order {M/TZY 
and (AIJTZ)^ displayed in Eq. (4.5 1 are not required in 
the determination of Q. 

we inte- 



Instead of proceeding directly with Eq. (4.9 



grate it over a cross-section of the horizon anc 
a differential equation for the surface area 



A{v) 



work out 



(4.11) 



Its rate of change is .4 = J ddd<j), and the second 
rate of change is ^ = J dvQ^d9dip+ J ^/^ dOdip; here 
we may once more ignore the second integral involving 
62. 

The differential equation that governs the behavior of 
the area function is 

KqA- A= [ aABCr'^^Vldedcjy + Oid). (4.12) 



The steps required in the evaluation of the right-hand 
side are described in Appendix O and the end result is 



kqA — A — SttJ-', 



(4.13) 



where 



16 



16 
4725' 



(c? cab 



16 



Sabc^'^"^ +0(9) 



(4.14) 



might be called the flux Junction. We notice that T is 
positive-definite, and that it cannot be expressed as a 
total derivative with respect to v. 



As we show in Appendix O the solution to Eq. (4.131 



IS 



""i(«)=^(.) + ^(") ' 



0(9). 



(4.15) 



We observe that the second and third terms on the right- 
hand side are total time derivatives, and that it is natural 
to bring them over to the left-hand side of the equation. 
We write our final result as 



Stt 



A* = J" + 0(9), 



(4.16) 



in terms of a "renormalized area function" defined by 
A* = A~%ttT / kI^BttJ^ / nl + OiW^^ /n^). The shift in 
area is of order /TZ^ and corresponds to a correction 
to the relation r = 2M of fractional order (M/TZ)^; this 
is well beyond the level of accuracy of Eq. (4.1). 



Equation (4.16) describes the rate at which the renor- 
malized black-hole area .4* increases as a result of the 
tidal interaction. We refer to this phenomenon as tidal 
heating, recalling the deep analogy that exists between 
the relativistic tidal dynamics of black holes and the New- 
tonian tidal dynamics of viscous bodies [Ml |7D] . The 
phrase also recalls the fact that in black-hole thermo- 
dynamics, the area plays the role of entropy S, while 
the surface gravity plays the role of (Hawking) temper- 
ature T; in this context we write {kq/8'k)A* ~ TS and 
interpret Eq. ( 4.16| as describing a flow of heat across 
the horizon. By the first law of black-hole mechanics, 
dM = {K/8TT)dA, the equation also describes the rate 
at which the black-hole mass increases during the tidal 
interaction. 



V. BACKGROUND METRIC: DERIVATION 

Our main results were presented in the preceding sec- 
tions, and we now turn to a detailed derivation of these 
results. In this section we provide a derivation of the 



3.3) and (3.4); an alterna- 



background metric of Eqs. 
five derivation is sketched in Appendix |D[ In Sec. VI 



present a derivation of the black-hole metric of Eqs. (3.5 1 
and (|3Jt. 



we 



A. Kinematical properties of the metric 

We construct the metric of a vacuum region of space- 
time that surrounds a timelike geodesic 7, and we adopt 
the light-cone coordinates {v,r,9^). These are centered 
on the world line, and the geometrical meaning of the co- 
ordinates was specified in Sec. |III A[ The three defining 
properties listed there give rise to four important condi- 
tions on the metric. 
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To spell them out we introduce the dual vector := 
—daV = [—1, 0, 0, 0], which is normal to hypersurfaces of 
constant advanced time v. According to the first prop- 
erty, this dual vector is null — g°'^Ca^i3 = — and this 
condition immediately implies that g""^ — 0. The vec- 
tor is tangent to the null generators of the light cones 

V = constant, and the second and third properties imply 
that its components in the light-cone coordinates must be 
given by r = [0, -1, 0, 0]; that T = €^ = means that 

V and are constant on the generators, and C = 
means that r is an affine parameter that decreases as the 
generators converge toward the world line. The relation 

= 5"^^/3 then implies that the inverse metric must 
satisfy the conditions 



ff™ = = 



5- = 1. 



(5.1) 



Calculating the inverse, we find that the metric must 
satisfy 



gv 



= 1. 



(5.2) 



These statements are exact, and follow from the light- 
cone nature of the coordinates. The other nonvanishing 
components of the metric are g^^, gvA, and gAB- We 
have not yet made a choice of normalization for the ra- 
dial coordinate r, nor a choice of axes for the angular 
coordinates 6^. 

From the quasi-spherical coordinates (r, 9^) we con- 
struct a system of quasi-Cartesian coordinates a;° = 
ril.°'{0^) and transform the spatial components of the in- 
verse metric according to the rules described in Sec. |IIE[ 
We find that g"" = fi", but the form of is not con- 
strained by the conditions of Eq. ( |5.1[ ) . 

At this stage we impose the additional condition that 
the inverse metric should be locally flat near the world 
line, so that g"'^ = 5°'^ + h"-^ , with h°'^ going to zero 
(as r^) when r — s- 0. This condition implicitly specifies 
a normalization for r, which reduces to the usual Eu- 
clidean distance in the immediate vicinity of the world 
line. And it specifies a set of reference axes for the an- 
gular coordinates, which are aligned with the Cartesian 
directions associated with the coordinates x°'\ the fact 
that scales as (instead of r) when r ^ implies 
that the world line is unaccelerated and that the axes are 
locally nonrotating. 

Our expression for the inverse metric in quasi- 
Cartesian coordinates is therefore 



7™ = 0, 



cab I lab 



(5.3a) 
(5.3b) 
(5.3c) 



with 



The last equation indicates that as r — )■ 0, 
^-2f^ab approaches the (direction-dependent) limit h"^^, 
r~^{r~^h°'^ — hf') approaches hf', and so on. 



To proceed it is useful to introduce a decomposition of 
h"^^ into longitudinal and transverse pieces. We write 

h''^^n''n^A + n''A'' + A''n'' + A'''' (5.5) 



with 



flaA" = = riaA'''' 



(5.6) 



Here and below, indices on fi" are lowered with the Eu- 
clidean metric Sab- The first term in Eq. (5.5 1 is the longi- 
tudinal piece of h""^, and A := h°'^il.a^b is its component 
in the direction of the unit radial vector fi". The last 
term is the transverse piece of /i"'', and := ^""^1^ Jt^'^'^ 
are the components in the directions orthogonal to fi"; 
7"^, is the projector of Eq. (2.151. The second and third 
terms are the longitudinal-transverse piece of /i*^^, and 
j^a ._ 'yO-^fi^b^^ contains the relevant components. The 
six independent components of h"^'' are contained in A 
(one component). A"- (two independent components), 
and A"^^ (three independent components). It is some- 
times useful to further decompose into trace and 
tracefree pieces, but we choose not to do so at this stage; 
we will find in due course that the Einstein field equations 
auto mati cally enforce SabA"'^ ~ 0. The decomposition of 
Eq. (5.5 1 can be ap plied individually to each h'^ that 
appears in Eq. (5.4 1; this defines An , Af , and A'}^ 



We write the inverse metric of Eq. ( 5.3 1 as g 



h"^ , with rj"^ denoting the inverse of the Minkowski met- 
ric in light-cone coordinates (with components ry"" = 0, 
7]"" = f]", and r]"'' = S"'') and h"'^ = 0{r^) denoting a 
perturbation (with components h'"'" — 0, h^'^ = 0, and 
h"'''). The metric is then g^p — rj^p — h^p + ha^h'^p + 
O(r^), where all indices are lowered with the Minkowski 
metric rjap (with components T^^y = — 1, rjya — ^a, and 
"Hab — lab)- A Straightforward calculation using Eq. (5.51 
reveals that 



gvv = -1 - A + AaA" + 0(r5), 
5™ = ^la-Aa + AabA^ + 0{r^), 
gab ^ lab - Aab + AacA\ + 0{r') . 



(5.7a) 
(5.7b) 
(5.7c) 



In these expressions it is understood that indic es o n A" 
and A°^ are lowered with 6ab\ by virtue of Eqs. (5.6 1 this 



operation is equivalent to lowering indices with riab — 
lab- It is also understood that A, and A"'^ can be 
expanded in powers of r as in Eq. ( |5.4[ ) ; we have 

A = r^A2 + r^A:i + r^Ai + 0{r^), (5.8a) 
Aa = r^A^a + r-3^3a + A^a + 0(r5), (5.8b) 
Aab = r^A2ab + r^A^ab + r^Aiab + O(r^), (5.8c) 

and a term like AaA^ reduces to r^A2a^2 + O(r^). 

We now return to the quasi-spherical coordinates 
(r, 6'"^). We define the angular version of the vector po- 
tential Aa by 



A, 



Aa^% 



(5.9) 
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^ .— dA^"', and we define the angular version of 



where ff^ 
the tensor potential Aab by 



(5.10) 



In spite of the suggestive notation, these are not the com- 
ponents of the Cartesian tensors Aa and Aab in spherical 
coordinates; for these we have Aadx°'/dr = Aa^"" — 
and Aadx"^ / d9^ = rAa^'A = "tAa, with similar equa- 
tions holding for Aab- 

Using the identities of Eqs. (2.18) we find that the 
spherical-coordinate form of the metric is 



9v 



9vA 



-l~A + AAA'^ + 0{r'^ 



1, 



-rAA + rAABA"" + Oir'^), 



(5.11a) 
(5.11b) 
(5.11c) 

gAB = r^^AB - t^Aab + t^AacA^ + O(r^). (5.1 Id) 

It is understood that in these expressions, indices on A a 
and Aab are raised with . As we saw previously, 
each potential can be expanded in powers of r, so that 

A = r^A2 + r^A:i + r'^Ai + 0{r'^), (5.12a) 
Aa = r'A^A + r^AsA + A^a + 0{r'), (5.12b) 
Aab = r'^A^AB + Aj,ab + t^A^ab + 0{r^). (5.12c) 

The coefficients A^ AnAi and AnAB are assumed to de- 
pend on 9"^ only, so that the dependence of the metric 
on r is contained explicitly in the power expansion. 



B. Field equations I 



The metric forms of Eqs. (5.7 1 and (5.11) follow di 



rectly from the light-cone nature of the coordinate sys- 
tems, and they embody the purely kinematical require- 
ments imposed by the choice of coordinates. To obtain 
more information we must impose the Einstein field equa- 
tions. 



In this first stage we involve the metric of Eqs. (5.11) 



in a computation of R"'" , the time-time component of 
the Ricci tensor. Setting this to zero order-by-order in r 
reveals that 



^^^AiAB ~ ^A2abA2^ 



(5.13a) 
(5.13b) 

(5.13c) 



n^''A2AB - 0, 

n^^'AsAB = 0, 

3 

5' 

In quasi-Cartesian coordinates these equations read 

7°'^2a6 - 0, (5.14a) 
l^'Asab = 0, (5.14b) 



7° ^4ab — -zA2abA2 

5 



ah 



(5.14c) 



These equations mean that Aab is tracefree through order 

\ab 

5 ' ^ ^zao^i2 



r'^, and that its trace is given by ^""^ Aab — \t'^ A2abA'^^ - 



We use this observation to simplify the form of the 
metric. Continuing to work in the quasi-Cartesian co- 
ordinates x", we re-express A/^^ab as a sum of trace and 
tracefree pieces: 



A 



4ab 



A 



iab 



3 

10 



labA2cdA2^, 



(5.15) 



where the new A^ab, like A2ab and A^ab, is now known to 
be tracefree. We also exploit the identity 



1 



A2acA2b = ^JabA2cdA2 , 



(5.16) 



which holds for any symmetric-tracefree tensor A2ab- 
Making these substitutions in Eq. (5.11), we arrive at 



-1 - r'A,, 



r^A3 - r^Ai 



r^A2aA-2 



9vv = -i - r ^2 

+ 0(r5), (5.17a) 

gva ^^a~ r^A2a - Asa - A^a + A2abA. 

+ 0(r5), (5.17b) 



qab = lab - r^A2ab ~ A^iab 



'A 



4a f) 



\r'labA2adAf + 0{r''). 



(5.17c) 



It is understood that the vector potentials A2a, ^3a, ^4a, 
and A2abA\ are transverse, in the sense that they are all 
orthogonal to fi". And it is now understood that the 
tensor potentials A2ab, A^ab, and A4^ab are &oi/i transverse 
and tracefree; the potential 7a6^2cd^2'' transverse and 
pure trace. 

In quasi-spherical coordinates the metric is 



: -1 - r^A2 - r^Aa - r^A^ + r^A2AA^ 
+ 0(r5), (5.18a) 
: 1, (5.18b) 
: ~r^A2A - r^A-iA ~ t^A^a + r^AiABA^ 
+ 0{r^), (5.18c) 

QAB = r'^^AB - r'^A2AB ~ r^A'iAB - r^AiAB 



gv 



9vA 



\T^nABA2CDA^'' + 0{r''). 



(5.18d) 



Here also the tensor potentials are tracefree, except for 
the term proportional to ^ab, which is pure trace. 



C. Field equations II 



The metric of Eqs. (5.17) and (5.18) is a partial solu 



tion to the Einstein field equations. In its quasi-Cartesian 
form the metric involves the transverse potentials A2a, 
A^a, Aia, ^2ab, ^3o6, and A4^ab, and imposing the vac- 
uum equation R""" = has revealed the important fact 
that the tensor potentials are all tracefree. In the quasi- 
spherical form of the metric, the potentials have compo- 
nents in the angular directions only, and these depend 
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on 9"^ only. To obtain a complete solution to the field 
equations we must now determine the potentials. 

We rely on Zhang's work [TB] , which allows us to state 
that the metric of a vacuum region of spacetime around 
a timelike geodesic is a functional of two (and only two) 
sets of tidal moments £aia2---ai o,nd Baia2---au tidal 
moments are STF tensors that depend on proper time 
on the world line, and they are related to components of 
the Weyl tensor (and its derivatives) evaluated on the 
world line. For a complete description of the metric one 
requires an infinite number tidal moments; for an ap- 
proximate description one requires a finite number. In 
our case the construction of the metric shall involve the 
quadrupole moments £ab and Bab, the octupole moments 
Sabc and Babci and the hexadecapole moments Sabcd and 
Babcd- These were introduced in Sec. II A and their scal- 



ing properties are described in Sec. II B 



Collecting the observations summarized in the preced- 
ing two paragraphs, we obtain three important guiding 
rules for the construction of the metric: 

1. The metric is constructed from scalar potentials, 
vector potentials that are transverse, and tensor 
potentials that are transverse and tracefree. 

2. The potentials depend on the angles 9"^ only; they 
are independent of r, which appears as a multiplica- 
tive factor in front of the potentials. 

3. The potentials depend on two sets of tidal mo- 
ments. 

The rules imply that the metric must be constructed 
from the tidal potentials introduced in Sees. 



II D and 



HE There are no other possible building blocks for the 
metric. 

We begin the construction of the metric with the de- 
termination of A2, A2A, and A2AB- (In practical mat- 
ters it is convenient to deal with the quasi-spherical rep- 
resentation of the metric, because the r-dependence is 
then explicitly known.) These terms occur at order 
in gyy, r^ in gyA, and r'* in gAB, and proper dimension- 
ality requires that the potentials scale as TZ~^. Equa- 
tion (2.121 then implies that the potentials must be con- 
structed from Sab and Bab- The possible building blocks 
are listed in Tables VIII and |IX| and we write A2 = a£ 
A2A 



b£\ + pB\, A 



12AB 



cE' 



AB 



qB 



AB^ 



where 



b, c, p, and q are undetermined numerical coefficients 
From this we form the metric g^ = — 1 — r^A2 + 0{r^) 



gvr — 1 
.4A 



gvA = -r"A2A + 0(r''), and gAB = r^^AB - 
r"A2AB + 0{r^), which we substitute into the vacuum 
field equations. (At this stage of the computations the 
D-dependence of the tidal potentials can be ignored, be- 
cause their w-derivatives are suppressed by a factor of 
order 7?.^^ relative to the spatial derivatives. Another 
simplifying move is momentarily to switch off the (fy 
dependence of the metric by adopting f g and B}^ as the 
only nonvanishing components of the tidal moments.) 
The exercise returns the relations h = |a, c = |a, and 
q = \p, which leaves a and p as undetermined constants. 



To obtain a and p we compute the frame components 
Caobo and CabcQ of the Weyl tensor in the limit r 0, 
and demand that the results agree with Eqs. (2.4 1. From 
this we find that a = 1 and p = — |, and the metric is 
known through order TZ^^ . 

For the computation of the frame components of the 
Weyl tensor we need 



- 1, 

na 



0, 



0, 



na 



(5.19a) 
(5.19b) 



the components of the tetrad vectors in the light-cone 
coordinates. Here fi;^ := ^atfi'^^ri^, and the factor of 
compensates for factors of r in the angular compo- 
nents of the metric. It is easy to show that the vectors 
are orthonormal and parallel-transported along in the 
limit r — !■ 0. We observe that while some components of 
the Weyl tensor are ambiguous or go to zero in the limit 
r — > 0, the frame components are well-behaved and have 
a well-defined limit. 

We continue with the determination of A3, A^a, and 
AzAB- Here the potentials must scale as TZ"^ , and 
Eq. 

Eab 1 



of 8\b, 



(2.121 implies that they must be constructed from 
Bab , Eabc , and Babe ■ The possible building blocks are 
listed in Tables|VIIl]and|lX] and we express A3 as a linear 
combination of and £°, A^a as a linear combination of 
and B\, and A^ab as a linear combination 

, £AB^ ^°AB- The numerical coefficients 
are determined in two steps. First, we construct an im- 
proved metric by appending the terms of order TZ~^ , and 
we substitute it into the vacuum field equations. (In this 
step the u-dependence of the quadrupole tidal potentials 
must be taken into account, again keeping in mind that 
the w-derivatives are suppressed by a factor of order TZ~^ 
relative to the spatial derivatives. It is still possible mo- 
mentarily to switch off the </)-dependence of the metric 
by adopting ^q, Bq, £q, and Bq as the only nonvanishing 
components of the tidal moments.) The first step leaves 
two coefficients undetermined, a common multiplicative 
factor in front of Sabc, and another common factor in 
front of Babe- These are determined in the second step, 
in which we compute the frame components C'aobQ\c and 
CabcO\d of the covariant derivatives of the Weyl tensor in 
the li mit r — s- 0, and demand that the results agree with 
Eqs. (2.6 1. The metric is now known through order TZ~^. 



We complete the computation of the metric with the 
determination of A4, A4A, and Aaab- Here the poten- 
tials must scale as TZr'^, and Eq. (2.121 implies that they 
must be constructed from £ab, Bab, £abc. Babe, £abcd, and 
Babcd- That is not all, however, because at order Tlr'^ 
we must also include potentials that are generated by 
quadratic combinations of Eab and Bab', the list of po- 
tentials is long, and this makes the computations much 
more involved than in the previous cases. The possible 
building^blocks are listed in Tables [Vml |IX] |X] |x!l [XH] 
We express A4 as a linear combination of 5"', 
V\ V\ Q", Q^ Q\ G', and G". We express 



and |XIII 
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A4A as a linear combination of £\, £\, £\, B\, B\, B\, 
V\, V\, Q\, Q\, g%, g%, n\, and n\. And we express 
^iAB as a linear combination of ^°ab^ ^ab^ ^ab^ 

^AB^ ^AB' T^AB' 'PaB' Q-AB' Q-AB' ^aB' ^AB^ ^AB> 

and 'H\b- addition to all this the angular compo- 
nents of the metric contain a bilinear term proportional 
to A2AB^^^ = which can 



solution 



be simplified with the help of Eqs. (B2|. The many nu- 



merical coefficients that appear in the metric at order 
TZ~^ are determined as we did previously. Substitution 
of the metric in the vacuum field equations determines 
all but two coefficients, the common multiplicative fac- 
tors in front of Eabcd and Babcd- These are determined by 



computing the frame components Cai 



ObO\cd 



and C, 



abcG\de 



of the second covariant derivatives of the Weyl tensor in 
the limit r 0, and demanding that the results agree 
with Eqs. (2.7 1. The metric is now known throug h or - 
der TZ~'^, and the end result is displayed in Eqs. (3.4 1. 



The quasi-Cartesian representation of Eqs. (3.3 1 can be 
obtained directly from this. 

This completes the derivation of the background met- 
ric. An alternative derivation is presented in Ap- 
pendix |D] 



ds'^ = -fdv'^ + 2dvdr 



■dn^ 



(6.1) 



where / : 

,2 n JJ,2 



1 - 2M/r and dn'^ = ^ABde^de^ = dO^ + 
sm't/CLfp". To construct the perturbation we rely heav- 
ily on the formalism of Martel and Poisson [33] , and we 
implement the light-cone gauge of Preston and Poisson 
[51] . The computations that lead to the black-hole metric 
are extremely lengthy, and are presented in four highly 
technical subsections. In Sees. IVI Al and IVIBI we con- 
struct the linear piece of the metric perturbation, and in 
Sees. I VI C| and |VI D| we turn to the bilinear piece. The 
results are collected in Sec. IVI El 



A. Linear perturbation: Preparation 

1. Form of the metric perturbation, and considerations of 
gauge 



VI. BLACK-HOLE METRIC: DERIVATION 



We next turn to a derivation of the black-hole metric of 
Eqs. (3.7 1, a vacuum perturbation of the Schwarzschild 



The metric of Eqs. (3.4 1 specifies the asymptotic condi- 
tions (when r ^ 2M) for the metric of a tidally deformed 
black hole. Focusing our attention on the linearly per- 
turbed piece of the metric, we write it as 



9vv = -/ - r^elf^ + - ^r^elS' - \r^^\^° + \r^(^\t - \^r^e\e + bilinear, (6.2a) 

gvr = 1, (6.2b) 
9.A = -\r\^,£\ - fe^S^) + \r\elt\ - 6^56^) - - blB\) - \r^el£:, - blB\) 

+ ^-r' {el£% - blB\) - ^r' {e\£), - b\B\) + bihnear, (6.2c) 

9AB = r'^AB - lr^{e^j£\B ^Wab) + " ^I^ab) " " blB^^) 

- Ir" {e^£\B - brBAB) + (e°8^AS - ^"s^as) - {e'l^AB - K^'ab) + bihnear, (6.2d) 



where the undetermined radial functions e5j(r), e^(r), 
e5j(r), 65j(r), 6^(7"), and (r) are all required to ap- 
proach unity when 2M /r — > 0. This ansatz for the met- 
ric is motivated by the facts that (i) it reduces to the 
Schwarzschild metric when the tidal fields are turned off; 
(ii) it reduces to the (linear piece of the) background 
metric when M — ^ (so that the radial functions all be- 
come equal to unity); and (iii) its expansion in terms of 
tidal potentials constitutes a decomposition of the metric 
perturbation into a complete basis of spherical-harmonic 
modes. An important aspect of the metric is that its v- 



dependence is assumed to be slow (with a time scale of 
the order of TZ) and contained in the tidal moments; we 
rule out time-dependent processes that take place over 
time scales comparable to 2M . 



We write the black-hole metric as gap — Qap + Pap, hi 
which gaj3 is the Schwarzschild metric of Eq. (6.1 1 and 
Pai3 is the tidal perturbation. In Eq. (6.2) the perturba- 



tion is presented in the Preston-Poisson light-cone gauge, 
which enforces the conditions p^r = Prr = PrA = 0. The 
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even-parity sector of the perturbation is 

_ ilmylm 
Fvv — / ^ "-vv ^ I 

Im 

PvA — 2_^Jv ^ A 1 

Ira 

PAB = Y.^K'^^abY''^ + G'^ri-g), 



(6.3a) 
(6.3b) 
(6.3c) 



Im 



where the reduced perturbations h}^, j^™, if'"*, and G'"* 
depend on v and r only. Preston and Poisson show that 
for vacuum perturbations, the gauge can be refined to 
also enforce if'™ = 0. This leaves /ij,™, j^™, and G'" 
as nonvanishing perturbations, and these can be read off 
from Eqs. (6.2). The odd-parity sector of the perturba- 
tion is 



PvA 
PAB 



/in y\ 



Im 



AB^ 



(6.4a) 
(6.4b) 



Im 



and /12™ depend on v and r only; these also 



in which h. 
can be read off from Eqs. (|6.2| 



Preston and Poisson show that the residual gauge free- 
dom that remains within the even-parity sector of the 
perturbation is a family characterized by an arbitrary 
function a''"(w) — one function for each / and m. Under 
such a gauge transformation the perturbations change 
according to 



M 

hvv hyy + l{l + l)^a — 2 
+ 2rd, 

jv^jv- + - f]a + 



2r2 



/(; + l) 



G^G~ -a 



r 1(1+1) 



a 

(6.5a) 
(6.5b) 

(6.5c) 



where the Im label was omitted for ease of notation. In 
addition, Preston and Poisson show that the residual 
gauge freedom that remains within the odd-parity sec- 
tor of the perturbation is a family characterized by an 
arbitrary function a'™(u) — one function for each I and 
m. Under such a gauge transformation the perturbations 
change according to 



r a, 



h2 - 2r^a. 



(6.6a) 
(6.6b) 



We shall use the residual gauge freedom to specialize 
the light-cone gauge to a horizon-locking gauge defined by 
the requirements Pw = Pvr — PvA = at r = 2M. Since 
the light-cone gauge already enforces pyr — everywhere, 
the horizon-locking gauge requires 



'Am 



■Im 

Jv 



: 



at r = 2M. 



(6.7) 



Because the residual gauge freedom is limited to two 
functions a'''^{v) and a'™(u), it may seem doubtful that 
the three conditions of Eq. (6.7 1 can be imposed. Nev- 



ertheless, we shall see that the Einstein field equations 
do allow the specialization of the light-cone gauge to the 
horizon-locking gauge. We saw in Sec. |IV A| that the 
conditions of Eq. (6.7 1 ensure that the horizon keeps its 
coordinate description r = 2M in the perturbed space- 
time. 

The horizon- locking gauge, together with the vacuum 
field equations, imply the existence of a useful constraint 
on the values of h^™ and /12" at r = 2M. We examine 
the perturbation equation = (in the notation of 
Martel and Poisson) and evaluate it at r = 2M. Making 
use of the statement hi™ — returns 



= - 



(;- l)(/ + 2) dh2 



SAP 



dv 



d dhy 
dv dr 



(6.8) 



Integration with respect to v, keeping r anchored at r 
2M, yields 



h'^{v,r = 2M) 



dk 



Im 



{l-l){l + 2) dr 



(6.9) 



r=2M 



The constant of integration was set equal to zero to re- 
spect the ansatz of Eqs. (6.2 1, which implies that there is 



no v-independent term in the perturbation. We shall re- 
fer to Eq. (6.9 1 as the horizon-locking constraint; it refers 



to the odd-parity sector of the perturbation only. 



2. Redefinition of tidal moments 

In addition to making use of the residual gauge free- 
dom, the form of the metric perturbation can be ad- 
justed by redefining the tidal moments Eab, £abc, Bab, 
and Babe according to £ ^ £ -\-piM£ -\-p2M'^£ and 
B ^ B-^ qiMB + gsM^B +■■■ . As we shall see below, 
these redefinitions have the effect of inducing changes in 
the radial functions ejj, e^, b'^, and b°^. 



3. Even-parity sector; I 



According to Eqs. (6.2 1 and Table VIII the even- 



parity, I — 2 piece of the perturbation is described by 



Jv ■ 
G 



-r'^eiE + ^-r^e2e - ^r^e:i£ + ■ ■ 

-\r^ei£ -t- \r^e^£ - ^r^e^E + 
3 6 63 

--r^erf -I- j^'^^egf - -r^e^E + 



(6.10a) 
(6.10b) 
(6.10c) 



in which e„ := e1^{r) and £ £'^{v). For ease of nota- 
tion we have omitted the label Im = 2m on the perturba- 
tion functions; the radial functions do not depend on m. 
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Under a residual gauge transformation generated by the 



function a{v) 



the radial functions change according to 



Ci ei + ci- 



62 



62 + 2ci- 



Af3 



63 ^ 63 



7 

2^^ 



,.3 
A/3 



3ei - 
r 

28 



4C2 



A/5 



A/4 



14C2 



A/5 



C3- 



M3 



A/4 



64 



64 - 61 - 



65 ^ 65 



3 r 



3 

A/4 

4 



ci- 



3 r 



7 A/3 1 



66 



67 



66 



A/3 



2 r 



A/5 

5 ' 

A/5 

5 



1 AT^ 



67-61- 



69 — > 69 + 



2 A/3 8 M4 
5 

28 



ci- 



762- 



27'V 



A/4 



0^3 

9 r 



A/5 

5~' 



(6.11a) 

(6.11b) 
A/6 

j,6 ' 

(6.11c) 
(6.11d) 

(6. lie) 

(6.11f) 

(6.11g) 

(6.11h) 

(6.11i) 



Here the residual gauge freedom was reduced from a 
functional family characterized by an arbitrary function 
a{v) to a three-parameter family (with parameters ei, 62, 
and 63). This loss of generality is a choice that is mo- 
tivated by the observation that the even-parity, I — 2 
piece of the perturbed metric should be driven by a sin- 
gle function £{v), so that a{v) should involve only £{v) 
and its derivatives. Below we shall seek choices for ei, 
62, and 63 that enforce the horizon-locking conditions 
61 = 62 = 63 = 64 = 65 = 66 = at r = 2A/. 

Implementing the redefinition £ — > 5 — |pii\/£ + 



i^P2M^£ + ■ ■ ■ in Eqs. (|6.10|) has the effect of chang- 



ing the identity of the radial functions 61 • • • 69. 
become 



They 



61 - 


61, 






(6.12a) 


62 - 


-> 62 H 


M 

Kpi — 61, 
r 




(6.12b) 


63 - 


^ 63 H 


7 A/ 

^ -^Pi — 62 4 
6 r 


A/2 

-P2^6i, 


(6.12c) 


64 - 


64, 






(6.12d) 


65 - 


65 H 


2 M 

3 r 




(6.12e) 


66 " 


66 H 


7 M 

^ ^Pi — 65 4 

8 r 


1 A/2 
2 


(6.12f) 


67 - 

68 - 


67, 

68 H 


2 M 
h -Pi — 67, 
5 r 




(6.12g) 
(6.12h) 


69 - 


69 H 


35 A/ 

h — Pl— 68 

54 r 


2 M2 

+ 0P2— 67. 
9 r"^ 


(6.12i) 



Below we shall seek choices for pi and p2 that enforce 
68 = 69 = at r = 2Ar; this does not alter the conditions 



already imposed by the horizon-locking gauge. With all 
these choices implemented, only 67 will be nonvanishing 
at r = 2A/. 



4- Even-parity sector; I = 3 



According to Eqs. (6.2 1 and Table VIII the even- 



parity, I — 3 piece of the perturbation is described by 



-^r3ei£-+ ^r462£ 
3 6 



G 



1 

12 
1 

18 



r464£ - 



r367£ - 



1 

18 
1 

20 



(6.13a) 
(6.13b) 
(6.13c) 



in which 6„ := e°„{r) and £ := £^{v). For ease of nota- 
tion we have omitted the label Im = 3m on the perturba- 
tion functions; the radial functions do not depend on m. 
Under a residual gauge transformation generated by the 
function a{v) = —-^ciM'^£ — ^62A/5£ -I- • • • , the radial 
functions change according to 



A/5 



61 



62 



64 



65 



61+61- 



62 



64 



65 



5 ^ A/4 
5 A/4 
A/5 



A/5 



ci- 



3 
5 

-62 F 

4 r5 
A/4 



3 r 



.5 

A/5 

5 



62- 



1 A/6 

^62^, 



67 ^ 67 - 61- 



68 



5 A/4 

4" 



27^^ r 



0^2 
9 r 



A/5 
5 ■ 



(6.14a) 
(6.14b) 
(6.14c) 
(6.14d) 
(6.14e) 
(6.14f) 



Here (as before) the residual gauge freedom was reduced 
from a functional family to a two-parameter family. We 
shall seek choices for 61 and 62 that enforce the horizon- 
locking conditions ei = 62 = 64 = 65 = at r = 2M. 

Implementing the redefinition £ —> £ — ^piM£ + ■ ■ ■ 
in Eqs. (6.131 has once more the effect of changing the 



identity of the radial functions. They become 



61 
62 
64 
65 
67 



61, 

62 
64, 

65 
67, 



M 
Pl — 61, 



3 M 
-Pi — 64, 

4 r 



5 A/ 
-Pl— 67. 
9 r 



(6.15a) 
(6.15b) 
(6.15c) 

(6.15d) 
(6.15e) 

(6.15f) 



Below we shall seek a choice for pi that enforces 63 = 
at r = 2A/; this does not alter the conditions already 
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imposed by the horizon-locking gauge. With all these 


67- 


67 H 


h fci 


M2 


(6.19d) 


2M. 














68- 








(6.19e) 


5. Even-parity sector; I = 4 


69- 


69 H 




„4 ■ 


(6.19f) 



According to Eqs. (6.2) and Table VIII the even- 
parity, I = 4 piece of the perturbation is described by 



hyy 

jv = 

G = 



1 4 p 

r eiC 

12 ^ 



1 

60 
1 

120 



(6.16a) 
(6.16b) 
(6.16c) 



in which e„ := e''j^{r) and £ := £^{v). For ease of nota- 
tion we have omitted the label Im = Am on the perturba- 
tion functions; the radial functions do not depend on to. 
Under a residual gauge transformation generated by the 
function a{v) = —-^ciM^S + • • • , the radial functions 
change according to 



M6 



ei 



64 



ei+ci- 



64 ■ 



67 - Ci 



9 Af5 

-ci — ^ 
4 r5 



1 



(6.17a) 
(6.17b) 
(6.17c) 



We shall seek a choice for ci that enforces the horizon- 
locking conditions ei = 64 = at r = 2M. This leaves 
67 as the only nonvanishing function at r = 2M. 



Odd-parity sector; I = 2 



According to Eqs. (6.2) and Table IX the odd-parity, 



I — 2 piece of the perturbation is described by 

4 



3 6 



63 



r%6B + 



h2 = lr%B-^r''bsB+l^r%B 



(6.18a) 
(6.18b) 



in which 6„ :— b1^{r) and B B'^{v). For ease of nota- 
tion we have omitted the label Im = 2m on the perturba- 
tion functions; the radial functions do not depend on to. 
Under a residual gauge transformation generated by the 
function a{v) = -^kiAPB-^ ^k^M^B- ^k3M^B+-- ■ , 
the radial functions change according to 



fci 



(6.19a) 
(6.19b) 

(6.19c) 



Here the residual gauge freedom was reduced from a func- 
tional family characterized by an arbitrary function a{v) 
to a three-parameter family (with parameters ki, k2, and 
fca). This loss of generality is a choice that is motivated 
by the observation that the odd-parity, / = 2 piece of 
the perturbed metric should be driven by a single func- 
tion B{v), so that a(v) should involve only B{v) and its 
derivatives. 

Below we shall seek choices for fci, fc2, and k^ that 
enforce the horizon-locking conditions 64 = 65 = fee = 
at r = 2M . Since 64 is gauge invariant, this will be 
achieved if and only if automatically vanishes at r = 



2M. Using Eqs. (6. 18 1 we find that the horizon-locking 
constraint of Eq. ( |6.9[ ) implies 

&7(2M) = -Mh'^{2M), 
3 

1 

5 
4 

'9^ 

in which a prime indicates differentiation with respect to 
r. The last equation will be used to determine k^. 

Implementing the redefinition B ^ B — ^qiAIB + 



&8(2M) = 
&9(2M) = 



-Af6^(2M), 
[m6^(2M), 



(6.20a) 
(6.20b) 

(6.20c) 



j^q2M^B + ■ ■ ■ in Eqs. (|6.18|) has the effect of chang- 



ing the identity of the radial functions 64 • 
become 



be 
67 

69 



64, 

b5- 



qi — 04, 



21 

be + T^gi — 05 
16 r 

br, 

08 + ^91 — 07, 
5 r 

35 M , 

36 r 



12-^ 



M2 



9 



'72- 



-67. 



. They 

(6.21a) 
(6.21b) 

(6.21c) 
(6.21d) 
(6.21e) 

(6.21f) 



Below we shall seek choices for qi and (72 that enforce 
&8 = ^9 = at r = 2M; this does not alter the conditions 
already imposed by the horizon-locking gauge. With all 
these choices implemented, only 67 will be nonvanishing 
at r = 2M. 



7. Odd-panty sector; I 



According to Eqs. (6.2) and Table IX the odd-parity, 
I — 3 piece of the perturbation is described by 



— r%B - —r'^bsB 
12 18 



(6.22a) 



27 



1 

18 



h2 = i-r^bjB - ^r^bsB 



1 

20' 



(6.22b) 



in which 6„ :— b°^{r) and B :— B°„^{v). For ease of nota- 
tion we have omitted the label Im — 3m on the pertur- 
bation functions; the radial functions do not depend on 
m. Under a residual gauge transformation generated by 

the function a{v) = -^hM^B + ^k2M^B H , the 

radial functions change according to 



64- 


bi, 




(6.23a) 


65- 




1, M3 

2 r'^ 


(6.23b) 


br- 


^bj 


+ ki-^, 


(6.23c) 


bs~ 


-> bs 


, M4 
+ fc2 — . 


(6.23d) 



Here the residual gauge freedom was reduced to a two- 
parameter family. 

Below we shall seek choices for fci and k2 that enforce 
the horizon-locking conditions 64 = 65 = at r = 2M. 
Since 64 is gauge invariant, this will be achieved if and 
only if 64 automatically vanishes at r = 2M. Using 



Eqs. (6.221 we find that the horizon- locking constraint 
of Eq. (6.9) imphes 



67(2M) = --Mb'i{2M), 





bs{2M) 



-^M6^(2M), 



(6.24a) 
(6.24b) 



in which a prime indicates differentiation with respect to 
r. The last equation will be used to determine fc2. 
Implementing the redefinition B ^ B — "^qiMB + • • • 



in Eqs. (6.221 has the effect of changing the identity of 



the radial functions. They become 



64 
bs 



&4, 
^8 



qi — 04, 



20 M 
27 r 



(6.25a) 

(6.25b) 
(6.25c) 

(6.25d) 



Below we shall seek a choice for qi that enforces bs — 
a,t r — 2M; this does not alter the conditions already 
imposed by the horizon-locking gauge. With all these 
choices implemented, only 67 will be nonvanishing at r = 
2M. 



8. Odd-parity sector; I = 4 



According to Eqs. (6.2) and Table IX the odd-parity, 
I — 4 piece of the perturbation is described by 



1 

60 



r%iB - 



(6.26a) 



h2 = r%B 

120 



(6.26b) 



in which 6„ := b)^{r) and B := B^„^{v). For ease of nota- 
tion we have omitted the label Im — 4m on the perturba- 
tion functions; the radial functions do not depend on m. 
Under a residual gauge transformation generated by the 
function a{v) = —■^kiAI'^B -I- • • ■ , the radial functions 
change according to 



&4, 



b7 



ki 



M4 



(6.27a) 
(6.27b) 



Here the residual gauge freedom was reduced to a one- 
parameter family. The horizon-locking condition is 64 — 
at r = 2M, and this will be achieved if and only if 
64 automatically vanishes dX r — 2M . Using Eqs. (6.26) 
we find that the horizon- locking constraint of Eq. (6.9) 
implies 



&7(2M) = --Mb'^{2M), 

y 



(6.28) 



in which a prime indicates differentiation with respect to 
r. This equation will be used to determine ki. 



B. Linear perturbation: Field equations 

1. Even-parity sector 

We describe in detail the method by which the pertur- 
bation equations are integrated in the quadrupole case 
I = 2; the same strategy is employed for the other mul- 
tipoles. We rely on the formalism of black-hole pertur- 
bation theory outlined in Martel and Poisson [33 . The 
method unfolds in a number of steps. 

In the first step the forms h^^ = —r'^ei{r)£{ev) -I- • • • , 

= -\r^ei{r)£{ev)-\ , and G = -\r'^ej{r)£{ev)^ 

are inserted within the even-parity field equations. Here e 
is a book-keeping parameter that reminds us that deriva- 
tives with respect to v are considered to be small; in 
this first step the field equations are expanded to order 
e°, and all terms in £ are neglected. In the notation of 
Martel and Poisson, the equation Q^"" = is automat- 
ically satisfied, and the equations Q^'' — 0, Q'"'" = 0, 
and — give rise to a system of three independent 
differential equations for the three unknown radial func- 
tions; the other field equations are related to these by 
the Bianchi identities. The general solution to the sys- 
tem of equations is easily obtained, and it depends on 
three integration constants. Imposing regularity at the 
event horizon determines one constant, and removes all 
terms proportional to log(r — 2M) from the radial func- 
tions. One of the two remaining integration constants 
is an overall multiplicative factor that is chosen so that 
when r ^ 2M, the functions ei, 64, and 67 all approach 
unity. The remaining constant is equivalent to the pa- 
rameter ci that appears in Eqs. (6.11 1; it characterizes 
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the residual gauge freedom that is still contained within 
the light-cone class of gauges. This last constant can 
be selected by simultaneously enforcing ei — 64 = Oat 
r = 2M; this is the horizon- locking condition of Eq. (6.7 1. 
The functions ei, 64, and ej are now completely deter- 
mined, and the constant ci has been chosen. 
In the second step the forms 

= -r^ei{r)£{ev) + ^r^e2{r)£{ev) + ■■■ , (6.29a) 
Y''eiir)£iev) + ^r^e5{r)£iev) + ■■■ , (6.29b) 
-^r^e7{r)£{ev) + es{r)£ {ev) + ■■■ (6.29c) 



Table IXIVl 



Jv 
G 



are inserted within the even-parity field equations. The 
functions ei, 64, and ey are already known from the first 
step, and the goal of the second step is to determine 62, 
65, and eg. The overdot still indicates differentiation with 
respect to v; the terms in £ are therefore of order e, and 
higher derivatives are of order and higher. In the sec- 
ond step the field equations are expanded to order e, and 
all terms in £ are neglected. As in the first step, the 
equations Q^^ — 0, Q^^ — 0, and Q"" — give rise to a 
system of three independent differential equations for the 
three unknown radial functions. The general solution is 
easily obtained, and it depends on three integration con- 
stants. As before, imposing regularity at the event hori- 
zon determines one of these, and removes all terms pro- 
portional to log(r— 2Af) from the radial functions. One of 
the two remaining constants is equivalent to the param- 
eter C2 that appears in Eqs. (6.11 1; it characterizes the 



residual gauge freedom that is still contained within the 
light-cone class of gauges. The other constant i s equ iva- 
lent to the parameter pi that appears in Eqs. (6.121; it 



corresponds to a redefinition of the quadrupole tidal mo- 
ment £. The constant C2 is selected by simultaneously 
enforcing 62 = 65 = at r = 2M ; thi s is once more the 
horizon-locking condition of Eq. (6.7). The constant pi 
is selected by also enforcing es = at r = 2M. The 
functions 62, 65, and eg are now completely determined, 
and the constants C2 and pi have been chosen. 

In the third and last step Eqs. (6.10 1 are inserted 
within the field equations, which are expanded to or- 
der e^. The general solution to the system of differential 
equations for 63, eg, and eg is easily obtained, and as al- 
ways it depends on three integration constants. The two 
that remain after imposing regularity at the event hori- 
zon are equivalent to the parameters C3 and p2 that ap- 
pear in Eqs. (6.11 1 and (6.12 1. The constant C3 is selected 
by simultaneously enforcing e3 = eg = at r = 2M. 
The constant p2 is selected by also enforcing eg = at 
r = 2M. The functions e3, eg, and eg are now completely 
determined, and the constants C3 and P2 have been cho- 
sen. 

The task of solving the linearized field equations in 
the even-parity, I = 2 sector is now completed, and we 
proceed in the same fashion for the I = 3 and I — 4 
sectors. The radial functions obtained here are listed in 



2. Odd-panty sector 



Essentially the same strategy is adopted for the odd- 
parity sector of the metric perturbation. We describe the 
details for the quadrupole {I = 2) sector. 

In the first step we insert the expressions hy — 
lr^b4r)B{ev) + ■■■ and /12 = lr%{r)B{ev) + ■■■ within 
the odd-parity field equations, which are expanded to 



)rder 



In the notation of Martel and Poisson, the 



= and — equations give rise to a system of 
two differential equations for the two unknowns and 
67. The general solution to this system depends on three 
constants of integration. The first is determined by de- 
manding regularity at r = 2M. The second is an overall 
multiplicative factor that is chosen so that 64 and 67 ap- 
proach unity when r ^ 2M. The third and final constant 
of integration is equivalent to the parameter ki that ap- 
pears in Eq. (6.19); it characterizes the residual gauge 



freedom that is still contained within the light-cone class 
of gauges. We find that &4 automatically vanishes at 
r = 2M, and we set the value of ki by imposing the 
first horizon- locking constraint of Eqs. (6.201. The func- 
tions 64 and &7 are now completely determined, and the 
constant ki has been chosen. 

In the second step we insert the expressions 

= lr%{r)B{ev) - lr%{r)B{ev) + ■■■ , (6.30a) 
3 6 

h2 = K%{r)B{ev) - ^r%{r)B{ev) + ■■■ (6.30b) 

within the odd-parity field equations, which are expanded 
to order e. As in the first step, the general solution to 
the system of differential equations for 65 and bs depends 
on three constants of integration. The first is determined 
by imposing regularity a,t r — 2M. The second constant 
is equivalent to the gauge parameter ^2 that appears in 
Eqs. (6.19), and the third is equivalent to the parame- 
ter qi that appears in Eqs. (6.21); this corresponds to 
a redefinition of the quadrupole tidal moment B. We 
find that 65 automatically vanishes at r = 2M, and we 
choose /c2 so that the second horizon-locking constraint 
of Eqs. (6.20) is satisfied. Finally, we choose qi so that 
bs also vanishes at r = 2M. The functions 65 and bg are 
now completely determined, and the constants ^2 and qi 
have been chosen. 



In the third and final step we insert Eqs. (6. 18 1 within 



the odd-parity field equations, which are now expanded 
to order . As before the general solution to the system 
of differential equations for 6g and bg depends on three 
constants of integration. The first is set by imposing 
regularity at r = 2M. The second is equivalent to the 
gauge parameter ks that appears in Eqs. (6.19), and the 



third is equivalent to the parameter 52 that appears in 
Eqs. ( 6.21[ ). We find that 6g automatically vanishes at 
r = 2M. We choose k^ so that the third of Eqs. (|6.20|) is 
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satisfied, and q2 so that 69 also vanishes at r = 2M. The 
functions be and bg are now completely determined, and 
the constants and ^2 have been chosen. 

The task of solving the linearized field equations in the 
odd-parity, I — 2 sector is now completed, and we proceed 
in the same fashion for the I = 3 and / — 4 sectors. The 
radial functions obtained here are listed in Table IXIVI 



have that b\, b^, bl, b}, bl, bl, b\, &°2, ^4, b\, and b\ all 



vanish at r = 2M. In addition, the freedom described in 
Sec. |VI A| to redefine the tidal moments was exploited to 



force 



-81 '^9; "=-8' ^8' ^9' 



and 6g 



also to vanish at r = 2M. 



With these conditions, the only radial functions that do 
not vanish at r = 2M are ey, 67, ey, 5y, 6y, and bj; their 
values at r = 2M are listed in the caption of Table |XIV| 



3. Light-cone gauge and horizon-locking condition 

The perturbation obtained in this section is presented 
in the light-cone gauge introduced in Sec. VIA The 
gauge is completely fixed: The residual gauge freedom 
that was initially left over was removed by imposing the 
horizon-locking conditions of Eqs. (6.7 1 and (6.9 1. In 



terms of the even-parity radial functions, this means that 
e\, 63, 63, 64, 65, Cg, e°i, e°2, e\, e%, e\, and e\ all vanish at 
r = 2M. In terms of the odd-parity radial functions, we 



C. Bilinear perturbation: Preparation 

1. Form of the metric perturbation, and considerations of 
gauge 

We next move on to the bilinear piece of the perturbed 
metric, which is generated by the linear quadrupole 
terms. The relevant pieces of the perturbed metric are 



-/ + linear + ^r^p^V"' + qlOT) + ^r^9\Q" + y\p\V^ + q\Q- 
1 



g^A = linear - ^r^92GA + ^r^h\n\ + ^r\plV\ + llqlQ\) 



75 



(6.31a) 
(6.31b) 

(6.31c) 



8 32 
gAB = r^^AB + linear + —r^nABiplV" + q-^QT) + i^^^^ab 9^0' 



"'3'^AB 



- ^^''^AB 9lQ° + r^9\Q\B + 'y'^AsW + + r\p\V\^ + q\Q\^) 



(6.31d) 



The bilinear metric perturbation is decomposed into mul- 
tipole moments, and with the help of Tables X] |XI[ |XII 



and XIII 



it can be expressed as in Eqs 
terms of scalar, vector, and tensor harmonics. The 



(6.3 



and (6.41 



metric of Eq. (6.31 1 includes all the tidal potentials that 
can contribute at bilinear order, including V^, Ti.\ 



Vab^ Sab' ^^"^ "^AB t^^t were missing from Eqs. (3.4 1 



All radial functions are required to approach unity as 



r ^ 00, except for p\, g\, p\, q\, h^, and 



/13, which are 



expected to approach zero as a power of 2M/r. 

The form of the metric perturbation can be altered by 
a residual gauge transformation that preserves the light- 
cone nature of the coordinate system. Because the gauge 
freedom was completely exhausted in the linear problem 
(by enforcing the horizon- locking conditions), the free- 
dom that is left over in the bilinear problem is purely 
bilinear — the coordinate shifts are necessarily propor- 



tional to 88, BB, or 8B. Under these circumstances the 
treatment of bilinear gauge transformations is identical 
to the linear treatment, and the residual gauge freedom 
is described by Preston and Poisson [SI] . In this context 
the residual gauge freedom is wider than in the linear 
problem: We can no longer enforce the tracefree condi- 
tion ri'^^pAB = 0, which was maintained in the linear 
problem. 

According to Preston and Poisson, the residual gauge 
freedom in the even-parity sector of the metric perturba- 
tion is described by 



jv -> jv + fa- 
K -> K + 2a + 



2 1 



3M\- ^; 

]h-2b 

r J 



2M 



r^c, 



b, (6.32a) 
(6.32b) 



-b + l{l + l)c, (6.32c) 
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G ^ G - -a~2c, 



(6.32d) 



the value of at r = 2M. A detailed examination of 



where a(w), h{v), and c{v) are the generators of the 
gauge transformation. The transformation is specific 
to each Im mode, and the complete gauge transforma- 
tion is obtained by summing over the relevant multipoles 
a = {0,1, 2, 3, 4}). 

The residual gauge freedom in the odd-parity sector of 
the metric perturbation is described by 



K 
h2 



riy — r a, 
h2 - 2r^a, 



(6.33a) 
(6.33b) 



where a{v) is the generator of the gauge transformation. 
This transformation also is specific to each Im mode, and 
the complete gauge transformation is obtained by sum- 
ming over the relevant multipoles {I — {2,4}). 

As we did in the linear problem, we shall exploit the 
residual gauge freedom to specialize the light-cone gauge 
to a horizon-locking gauge defined by the requirements 



■Im 
Jv 







at r = 2M. 



(6.34) 



These conditions, together with the field equations, give 
rise to the same horizon-locking constraint as in the linear 
problem: 



hi^{v,r = 2M) 



8M^ 



dh 



Im 



{l-l){l + 2) dr 



(6.35) 



r=2M 



This comes about for the following reason: In the nota- 
tion of Martel and Poisson ^33j , we find that while P[j^ 
is no longer identically zero in the bilinear problem (be- 
cause the linear, quadrupole piece of the metric perturba- 
tion produces an effective energy-momentum tensor for 
the bilinear metric perturbation), its value turns out to 
be zero at r = 2M; the derivation leading to Eq. (6.35) 



is therefore the same as the one leading to Eq. (6.9) in 
the linear problem. 

In addition to this first set of horizon-locking con- 
straints, another set arises as a consequence of the even- 
parity perturbation equations, in particular the Ql^ 
equation (again in the notation of Martel and Poisson). 
When evaluated at r = 2M, this equation implies 

a„if'"(f , r = 2M) ^ 4MQ^^{v, r = 2M) (6.36) 

when we also impose the horizon-locking conditions 



T Im 



J 



lin. 



at 



2M. The right-hand side of this 



equation is not zero, and the bilinear (5[^ can be com- 
puted from the linear pieces of the metric perturbation. 
Calculation reveals that the contributions to (5[^'j that 
originate from terms involving £ab and Bab in the linear 
perturbation vanish at r = 2M , so that Q\l^{r = 2M) 
comes entirely from terms that involve £ab and Bab- The 
end result is that Q\l^{r = 2M) is equal to a quantity 
quadratic in £ab and Bab that is differentiated with re- 
spect to V. Equation ( |6.36 ) can therefore be integrated 
with respect to v, and this gives rise to conditions on 



Eq. (|6.36|) shows that these are equivalent to the horizon- 

(6.37) 



locking constraints 



5 

16 



and 



5 

16 



93 = 93 = 

on the value of the radial functions at r = 2M. 



(6.38) 



2. Redefinitions 

In addition to making use of the residual gauge free- 
dom, the form of the metric perturbation can be adjusted 
by redefining the black-hole mass parameter M and the 
tidal moments Sat, Bab, £abc, and Babe- We consider the 
changes 

M ^miM^£pg£P'i + ^m2M^Bp,BP'' , 

15 15 

(6.39a) 

2 2 

Sab £ab - l^m:iM'^£pi^a£\) - ^m4M'^Bp(^al3^^y 

(6.39b) 



Bab Bab + -m^M £p(^aB\y 

£abc £abc - '2mQMepq(^a£''b^''c) ' 
Babe ^ Babe 

where the parameters m„ 
and where the factors of 4r 

io 



(6.39d) 
(6.39e) 



are dimensionless constants, 
I , I , and 2 were inserted for 
convenience. The redefinition of M formally introduces a 
time dependence in the black-hole mass parameter; this 
can be ignored at the level of accuracy maint ained in this 
work, because according to the first of Eqs. (6.39), M — 
0{M^ /TZ^). The impact of the redefinitions on the radial 
functions will be examined below. We observe that there 
is no need to consider a change such as £abcd £abcd + 
mg£(^ab£cd) + 'm-iQB{abBcd) , because this redefinition does 
not involve M; this ambiguity can be resolved at the level 
of the background metric. 



The redefinitions of Eqs. (6.39) are compatible with 



the parity rules spelled out in Sec. II C The rules forbid, 
for example, the presence of a £B term in the shift in 
£ab, and the presence of £^ and B^ terms in the shift in 
Bab- Similarly, a shift in Babe is ruled out, because the 
combinations of £ and B that could be involved would 
violate the parity rules. 



3. Even-parity sector; I = 



It is easy to pick out the monopole piece of the bilin- 



ear metric perturbation from Eqs. (6.31 1; it involves the 
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tidal moments V"" := EpqE^'^ and Q" : = gpp gP'?. Examin- 
ing the gauge transformation of Eqs. (|6.32l, we find that 



the changes for I = concern h^^ and K only, and that 
those are generated by two functions of time, a,(v) and 
h(v). Restricting the gauge freedom as we have done in 
the linear problem, we write a — j^M'^{ciP"' +(1x0^) and 
h = Y^M^{c2V"'+d2Q"'), with c„, d„ denoting dimension- 
less constants. Neglecting time derivatives, substitution 
within Eqs. (6.321 and comparison with Eqs. (6.311 re- 



veals that the gauge transformation produces the follow- 
ing changes in the radial functions: 



(6.40a) 
(6.40b) 
(6.40c) 
(6.40d) 











Pi ~ 


-^Pi - 


15 


DCl^ +2C2— ^, 

15 Af5 


Pa - 


-^Ps- 


4 

^ , 


4 r'' 


Qi ' 


^ li - 


h2di— - 
15 , 


6^17^ + 2^2^ 
15 AfS 


93 - 


<l3 ' 


4 





The redefinition of M in Eqs. (6.39) produces an addi- 
tional change in both pi and gj: 



Pi ^ Pi + 2mi- 



91 ^ 91 + 2?Tl2 - 



„5 " 



(6.41a) 
(6.41b) 



Below we shall use the horizon-locking condition pi^O 
at r = 2M to determine ci , leaving p"^ and dependent 
on C2 and toi. We shall next impose the horizon- locking 
constraint = 5/16 at r = 2M to determine mi. These 
conditions leave C2 undetermined, and this is chosen so as 
to simplify the form of the radial functions. Similarly, we 
shall use the horizon-locking condition — at r = 2M 
to determine di, leaving and dependent on d2 and 
TO2 ■ We shall next impose the horizon- locking constraint 
q^ = 5/16 at r = 2M to determine 7712- These conditions 
leave d2 undetermined, and this is chosen so as to simplify 
the form of the radial functions. 



4- Even-parity sector; I = 1 

The dipole piece of the bilinear perturbation involves 
the tidal potentials G'' and CJ^ defined in Table XII 



The gauge transformation for I = 1 concern hy^, jy, 
and K only. It is generated by three functions of time. 



which we write as a 



15 



ciApg", b 



15 



C2A'Pg', and 



c = j^c^M'^Q'' , with Cn denoting dimensionless constants. 
(There is an abuse of notation here. The functions a, b, 
and c are specific to each mode ^ = 1, m = {0, Ic, Is}; the 
notation G'' therefore refers to each harmonic component 
of the bilinear tidal moment, as listed in Table XII ) Ne- 



glecting time derivatives, substitution within Eqs. (6.32) 



and comparison with Eqs. (6.31) reveals that the gauge 



transformation produces the following changes in the ra- 
dial functions: 



g{ ^ g{ + 2c2 



M6 



92 ^ 92 



93 ^ 93 



C2 



15 Af4 



A/5 
15 , 

+ ^(ci - C2) 



2 ^ 



Af6 

r-6 '■ 



Af5 



(6.42a) 
(6.42b) 
(6.42c) 



Below we shall use the horizon-locking conditions g{ — 
52 = at r = 2Af to determine C2, leaving dependent 
on ci and C3. We shall next impose the horizon- locking 
constraint 53 = —5/16 at r = 2M to determine ci. These 
conditions leave C3 undetermined, and this is chosen so 
as to simplify the form of the radial functions. 



5. Even-parity sector; I = 2 

The quadrupole piece of the bilinear perturbation in- 
volves the tidal potentials T"^ and Q^, as well as their vec- 
torial and tensorial counterparts; these are defined in Ta- 
bles |X] and |XI[ The residual gauge transformation is gen- 
erated by a = ^M^iciV-^diQ"), b = fA/5(c27"' -1-^2 2'), 
and c = ^M'^ic^V-i-daQ"), in which V and Q" stand for 
the harmonic components of the tidal moments listed in 
Tables |X] and |XI[ Neglecting time derivatives, substitu- 
tion within Eqs. (6.32) and comparison with Eqs. (6.31) 



reveals that the gauge transformation produces the fol- 
lowing changes in the radial functions: 



Pi Pl+ 2C2 

pI ^ pI + 15(ci 
„ 45 



, A/5 Af 6 

C2j— 30ci- 



P l P 3 



Pi Pi 



91 91 



3 r 



2d, 



Af4 

a' 

Af4 

-^^ 

Af6 

^6 ' 



15 



(3ci - C2) 



A/5 



92 92 



'-"{di 



d2) 



91 



, 45 M4 



91^91 



^d. 



A/4 



30 A/"6 
A/5 

— (3dl -<i2)^: 

4 r° 
A/5 



A/5 
15 



(6.43a) 
(6.43b) 
(6.43c) 
(6.43d) 
(6.43e) 
(6.43f) 
(6.43g) 
(6.43h) 



The redefinition of Eq. (6.39) produces the additional 
changes 



(6.44a) 
(6.44b) 
^ (6.44c) 



p\~ 


-p\ 


^ f2 


Pl~ 


-^Pl 


+ 5to3^/, 


Pl^ 


-Pi 


4 A/2 / 



32 



91 <l\ +™4^/ 



5 

11 



"^4 — /, 



ql^ ql+ ^"^4 

9 r 



r 

2^ 



21VP 



(6.44d) 
(6.44e) 
(6.44f) 



on ci, C3, and mg. We shall next impose the horizon- 
locking constraint 33 — 5/16 at r = 2M to determine 
Ci, and the additional condition g\ — Q aX r — 2M to 
determine mg. These conditions leave ca undetermined, 
and this is chosen so as to simplify the form of the radial 
functions. 



Below we shall use the horizon-locking conditions p\ — 
P2 — ai r — 2M to determine C2, leaving p^^ dependent 
on ci, C3, and m^. We shall next impose the horizon- 
locking constraint = 5/16 at r = 2M to determine 
ci and the additional condition P4 = at r = 2M to 
determine TO3. These conditions leave C3 undetermined, 
and this is chosen so as to simplify the form of the ra- 
dial functions. Similarly, we shall use the horizon- locking 
conditions ql = = at r — 2M to determine ^2, leav- 
ing q"^ dependent on di, d^, and 7714. We shall next im- 
pose the horizon-locking constraint g3 = 5/16atr = 2M 
to determine di and the additional condition ^4 = at 
r = 2M to determine 1714. These conditions leave ^3 un- 
determined, and this is chosen so as to simplify the form 
of the radial functions. 



7. Even-parity sector; I = 4 

The hexadecapole piece of the bilinear perturbation 
involves the tidal potentials V*' and Q'', as well as their 
vectorial and tensorial counterparts; these are defined 
in Tables |X] and |XI[ The residual gauge transforma- 
tion is generated by a = -iM^(ci7"' + diQ'), b = 
-\M^{c2V' + dsQ"), and c = -\M'^{c:iV'' + dsQ-), in 
which 7"' and denote the harmonic components of 
the tidal moments listed in Tables |X] and XI Neglect- 
ing time derivatives, substitution within Eqs. (6.32) and 



comparison with Eqs. (6.31 1 reveals that the gauge trans- 
formation produces the following changes in the radial 
functions: 



6. Even-parity sector; I = 3 

The octupole piece of the bilinear perturbation involves 
the tidal potentials Q°, Q\, and Q\g, which are defined 
in Table |XIlj The residual gauge transformation is gener- 
ated by a = jciM^G", b = |c2M^0°, and c = IcsM^G", 
in which Q" stands for the harmonic components of the 



tidal moments listed in Table XII Neglecting time deriva- 
tives, substitution within Eqs. (6.321 and comparison 



with Eqs. (6.31) reveals that the gauge transformation 



produces the following changes in the radial functions: 



9°i - 


9°i 


+ 2C2^, 


(6.45a) 


92 ~ 


9°2 


+ 5(ci-C2)^-10ci — , 


(6.45b) 


53 - 


-^9°3 


-45c3^-y(6ci-C2)^, 


(6.45c) 


9l ' 


^ 9l 


-4c3^-4ci^. 


(6.45d) 



The redefinition of Eqs. (6.39) produces the additional 
changes 



M , 

5i 5i + ™6 — I 1 
r 



92 92 + |"l6 ^ 



5i-> 54 +3-6-1/+^ 



4Af3 



(6.46a) 
(6.46b) 
(6.46c) 



Below we shall use the horizon-locking conditions g\ = 
52 = at r = 2M to determine C2, leaving dependent 





-p\ 


+ 2C2^, 




(6.47a) 


P2- 


-p'2 


+ 10(ci-C2)^-20ci 


A/6 


(6.47b) 




-P3 


- 150c3— p - 15(10ci - 


, A/5 


(6.47c) 




-P\ 


Af4 Af5 
+ 4c3^+4ci-^, 




(6.47d) 


9i- 




+ 2d2^, 




(6.47e) 


12- 


-q'2 


+ 10(di-d2)^-20d] 


A/6 


(6.47f) 




-53 


A/4 

- 150d3— r - 15(10di - 




(6.47g) 




-5^ 


A/5 

+ 4^3-^ +4di-^. 




(6.47h) 



Below we shall use the horizon-locking conditions p'l — 
P2 — at r — 2A/ to determine C2, leaving depen- 
dent on ci and C3. We shall next impose the horizon- 
locking constraint P3 — 5/16 at r = 2M to determine ci. 
These conditions leave C3 undetermined. We shall find 
that P4 ^ at r = 2M, but that its value is independent 
of C3 , which is chosen so as to simplify the form of the ra- 
dial functions. Similarly, we shall use the horizon- locking 
conditions q'l = q2 = at r = 2M to determine d2, leav- 
ing q^^ dependent on di and ^3. We shall next impose 
the horizon-locking constraint q^ — 5/16 at r — 2 A/ to 
determine di. These conditions leave c?3 undetermined. 
We shall find that q\y^ at r = 2M, but that its value is 
independent of ^3 , which is chosen so as to simplify the 
form of the radial functions. 
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Odd-parity sector; I = 2 



D. Bilinear perturbation: Field equations 



The quadrupole piece of the bilinear perturbation in- 
volves the tidal potentials T-L\ and 'H\b defined in Ta- 
ble |XIII[ The residual gauge transformation is generated 
by a = |fciM'*7i'', in which Ti,'^ stands for the harmonic 
components listed in Table XIII Neglecting time deriva- 
tives, substitution within Eqs. (6.331 and comparison 



with Eqs. (6.31) reveals that the gauge transformation 



produces the following changes in the radial functions: 



M4 



(6.48a) 
(6.48b) 



The redefinitions of Eq. (6.391 produces the additional 
changes 



K ^ K 



4 AP , 



(6.49a) 
(6.49b) 



The horizon- locking constraint of Eq. (6.351 gives rise to 



him) 



9 dr 



(6.50) 



r=2M 



We shall verify that the horizon-locking condition 



/i2(2M) = is automatically satisfied, and use Eq. (6.50) 



to determine fci. 
hl{2M) = 0. 



We shall then choose TO5 so that 



Odd-parity sector; I = 4 



The hexadecapole piece of the bilinear perturbation 
involves the tidal potentials Ti\ and Ti\g defined in Ta- 
ble |XIII The residual gauge transformation is generated 



by a = ^kiM^H'', in which H'' denotes the harmonic 
components listed in Table |XIII[ Neglecting time deriva- 
tives, substitution within Eqs. (6.33) and comparison 



with Eqs. (6.31) reveals that the gauge transformation 



produces the following changes in the radial functions: 



2fci 



(6.51a) 
(6.51b) 



The horizon- locking constraint of Eq. (6.35) gives rise to 



hl{2M) 



1 . .dh'o 



-M- 



dr 



(6.52) 



r=2M 



We shall verify that h2{2M) = 0, and use Eq. ( |6.52[ ) to 
determine ki. We shall find that does not vanish at 
r = 2M, but that its value is gauge invariant. 



1. General strategy 

The solution to the bilinear problem requires a general- 
ization of the perturbation formalism employed in the lin- 
ear problem. The perturbed metric takes the schematic 
form g = g + epi + €^P2 + 0{e^), in which g stands for 
the Schwarzschild metric, epi is the linear perturbation 
calculated previously, and e^p2 the bilinear perturbation 
that we now wish to obtain; e ~ TZ~^ is a perturba- 
tive book-keeping parameter. While epi is linear in the 
tidal moments £ab and Bab, £^P2 involves terms of the 
schematic form ££, £B, and BB; those are decomposed 
into multipole moments and expressed in terms of the 
potentials V",V\V\ Q", Q^ Q\ Q\ Q\n\ and W. 
Notice that the bilinear problem involves the quadrupole 
tidal moments £ab and Bab only; at order ~ 7?."^ there 
is no need to involve the octupole and hexadecapole mo- 
ments that also appear in the linear perturbation. 

The Einstein tensor for the perturbed metric takes the 



schematic form G[g\ -\- eGi[pi 



'Gi[p2 



'G2[P1 



0{e^), where G[g] = is the Einstein tensor of the 
Schwarzschild metric, and eGi[pi] = is the first-order 
perturbation created by the linear perturbation epi . The 
remaining terms appear at second order. The first con- 
tribution, e^Gi[p2], is the Einstein tensor generated en- 
tirely from the second-order perturbation e^p2 ; Gi is the 
same linear differential operator that was encountered in 
the linear problem. The second contribution, e^G2[pi], 
is generated by the first-order perturbation, and it origi- 
nates in the nonlinearities of the vacuum field equations. 

The bilinear perturbation problem consists of finding 
solutions to the field equations 



Gl[p2 



(6.53) 



These have the same formal structure as the linear field 
equations Gi[pi] = 0, except for the fact that there is 
a source term on the right-hand side. The problem is 
tractable because the differential operator on the left- 
hand side is the same as in the linear problem, and be- 
cause the source term on the right-hand side can be com- 
puted directly from the known solution epi. Notice than 
time derivatives of £ab and Bab can be ignored when com- 
puting Gi[p2] and G2[pi]; these contribute at order TZ~^ 
and higher, and they do not affect the field equations at 
order ~ 7^-4. 

The strategy to solve the bilinear field equations is the 
same as in the linear problem. First p2 is decomposed 
in scalar, vector, and tensor harmonics, a task that was 
already accomplished in Eqs. (6.31). Second, the linear 



differential operator Gi[p2] is allowed to act on the per- 
turbation, and the result is again decomposed in spherical 
harmonics. Third, the effective source term — G2[pi] is 
computed and decomposed in spherical harmonics. And 
fourth, the perturbation equations are integrated, one 
mode at a time. In practice we rely on the Martel-Poisson 
perturbation formalism [33| to implement this strategy. 
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2. Even-parity sector 

The even-parity sector of the perturbation p2 involves 
the terms in P"' , V , V , Q-" , , Q\ , and G"; these are 
decomposed as in Eqs. (6.3 1. The source terms — G'2[pi] 



are computed from pi (which, we recall, involves the po- 
tentials £'' and B"* only), and these also are decomposed 
as in Eqs. ( |6.3[ ). In the notation of Martel and Poisson, 
the reduced source functions are denoted Q"", Q^^ , Q^^ , 
Q'',Q'',Q\ and Q», in which we suppress the use of the 
Im labels. Examination of p2 and calculation of —G2[pi\ 
reveal that the relevant multipoles are / = {0, 1, 2, 3, 4}. 

The perturbation equations for I = decouple into two 
sets, the first proportional to V"" and involving the radial 
functions p1 and pg, the second proportional to Q" and 
involving the radial functions q" and gj. For each set 
of equations the general solution depends on three in- 
tegration constants; these were denoted {ci,C2,rni} and 
{cii,d2,TO2} in Eqs. (6.40) and (6.41). The constants c„ 
and dn are chosen so as to enforce the horizon-locking 
condition p"^ = q"^ = Q aX r = 2M. The constants mi 
and TO2 are chosen so as to simplify the form of the radial 
functions; we use this freedom to remove terms propor- 
tional to in Pi/f and qi/ f, and terms proportional 
to x~'^ in P3// and (73//. The end results are listed in 



Table ISS 

The perturbation equations for I = 1 involve the dipole 
potentials G'' and the three radial functions gjj. The gen- 
eral solution to the set of coupled differential equations 
for the radial functions depends on three integration con- 
stants; these were denoted {ci, C2, C3} in Eqs. (6.42 1. The 
constants ci and C2 can be chosen so as to enforce the 
horizon-locking condition = at r = 2M. This leaves 
C3 undetermined, and it can be chosen so as to simplify 
the form of the radial functions; we use this freedom to 
remove a term proportional to in g\l J ■ The end 
results are displayed in Table |XV[ 

The perturbation equations for / — 2 decouple into 
two sets, the first proportional to and involving the 
four radial functions p^, the second proportional to Q*" 
and involving the four radial functions q\. For each set of 
differential equations the general solution depends on five 
integration constants; one can be eliminated by remov- 
ing all terms proportional to log(r — 2Af) from the radial 
functions, and the others correspond to {ci, C2, C3, m3} 
and {di,(i2,d3,TO4} in Eqs. ( |6.43| and ( |6.44[ ). The con- 
stants ci, C2, and TO3 are chosen so as to enforce the 
horizon-locking condition = at r = 2Af, while di, 
(i2, and are chosen to enforce q\ — ^ air — 2M . This 
leaves C3 and undetermined, and those are chosen so as 
to simplify the form of the radial functions; we use this 
freedom to remove terms proportional to in both 
p\/ f and q^/ f ■ The end results are listed in Table XV 



The perturbation equations for I — 3 involve the oc- 
tupole potentials Q° and the four radial functions g"^. The 
general solution to the set of coupled differential equa- 
tions for the radial functions depends on five integration 
constants; one can be eliminated by removing all terms 



proportional to log(r— 2M) from the radial functions, and 
the remaining constants were denoted {ci, C2, C3, mg} in 
Eqs. (6.45) and (6.46). The constants Ci, C2, and mg can 



be chosen so as to enforce the horizon-locking condition 
g°j^ — Q at r — 2M. This leaves C3 undetermined, and 
it can be chosen so as to simplify the form of the radial 
functions; we use this freedom to remove a term propor- 
tional to x~^ in 52//- The end results are displayed in 



Table IXVl 

The perturbation equations for I = 4 decouple into 
two sets, the first proportional to 7"' and involving the 
four radial functions pj^, the second proportional to Q'' 
and involving the four radial functions q*^. For each set 
of differential equations the general solution depends on 
five integration constants; one can be eliminated by re- 
moving all terms proportional to log(r — 2M) from the 
radial functions, and another must be set so that the ra- 
dial functions (except for p'l and q\) approach unity as 
r —f 00. The remaining three are denoted {ci, 02,03} 
and {c?i, (i2, (is} in Eqs. (6.47). The constants ci and C2 
are chosen so as to enforce the horizon-locking condition 
Pi=P2=P3 = 0sttr — 2M, and we find that the value 
of p\ at r = 2M is independent of C3 (and therefore 
gauge invariant); we choose C3 so that a term propor- 
tional to x^^ is removed from P2/ f- Similarly, di and ^2 
are chosen so as to enforce the horizon-locking condition 
q'l = q2 = q'^ = at r = 2M, and we find that the value 
of q\ at r = 2M is independent of ds (and therefore gauge 
invariant); we choose ^3 so that a term proportional to 
x^^ is removed from (72//- The end results are listed in 
Table |XVl 



3. Odd-panty sector 

The odd-parity sector of the perturbation p2 involves 
the terms in H'' and H'; these are decomposed as in 
Eq. (6.4). The source terms — G2[pi] are computed from 



Pi and also decomposed as in Eq. (6.4). In the notation 
of Martel and Poisson, the reduced source functions are 
denoted P", P'', and P, in which we suppress the use 
of the Im labels. Examination of p2 and calculation of 
— G2[pi] reveal that the relevant multipoles are I = {2, 4}. 

The perturbation equations for I = 2 involve the poten- 
tials W and the radial functions /ij ^^'^ ^3- The general 
solution to the set of coupled differential equations for the 
radial functions depends on three integration constants; 
one can be eliminated by removing all terms proportional 
to log(r — 2M) from the radial functions, and the others 
correspond to {fci, m^} in Eqs. (6.48 1 and (6.49 1. We ob- 
serve that 



at 



2M, and we use Eq. (6.50) to 



set the value of ki. Finally, we choose 7715 to enforce the 
additional condition that h''^ ~ at r — 2M. The end 



results are listed in Table IXVl 

The perturbation equations for I = 4 involve the poten- 
tials T-C" and the radial functions /12 and /13. The general 
solution to the set of coupled differential equations for the 
radial functions depends on three integration constants. 
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One of these can be eliminated by removing all terms pro- 
portional to log(r — 2M) from the radial functions, and 
another must be set so that both /ij and /I3 approach 
zero as r —^ 00. The remaining constant is denoted fci in 



Eq. ( |6.51[). W e observe that /I'j = at r = 2M, and we 
use Eq. (6.501 to set the value of fci. The end results are 
hsted in Table EVl 



E. Conclusion 



The linear piece of the metric perturbation was com- 
puted in Sees . | VI ^ and [VI B | an d the bilinear piece was 
computed in Sees. | VI Cj and VI D[ The calculation is com- 
plete, and after collecting results we obtain the black-hole 
metric of Eqs. (3.7 1. The quasi-Cartesian representation 



of the metric can immediately be constructed from this, 
and the expressions are given in Eqs. (3.5 1. 
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APPENDIX A: DECOMPOSITION OF TIDAL 
POTENTIALS IN SPHERICAL HARMONICS 

We wish to express the angular version of the tidal 
potentials hsted in Tables [U |ll] [III] |IV] [V) and |VI 



pansions in scalar, vector, and tensor harmonics. Before 
we proceed we record the useful identities 



(Ala) 
(Alb) 
(Ale) 



These quantities were all introduced in Sec. |IIE| of the 
main text. 

The general structure of the even-parity tidal poten- 
tials is 

(A2a) 

^^:A,u....kS^^---Q!'\ (A2b) 

^ - ^i^i.^Acdk.-k^n^' ■■■^^'+ la,A^'\ (A2c) 

in which Akik2---ki is a constant STF tensor of rank I. It 
is not difficult to show that these satisfy the eigenvalue 
equations 

r'^-i'"^DcDdA^^^ + l{l + l)A^^^ ^ 0, (A3a) 
T^^^-'D.DaA^^ + + \)- 1]^« = 0, (A3b) 



0, 



(A3c) 



where Da is a projected differential operator that acts 
as follows on arbitrary tensor fields: DaT/j-^i,^... :— 

la Ibi Ib2 ^p-^qiq2 - - 

The transformed potentials are 



^A 

1 (0 _ 00a nb 



a'a' ^n'XAak2-k,n'''---^'" 



Aab ~ '^^A^B^abks - ki^ ^ • • • f2 ' 



(A4a) 
(A4b) 

-UabA^'K (A4c) 



and we wish to expand them in the even-parity harmonics 
of Eqs. (2.20). These satisfy the eigenvalue equations 



n'^^DcDoY^"' + l{l + l)Y'"' = 0, (A5a) 
n^'^DcDoYjl'' 1] yj" = 0, (A5b) 

n^^DcDoYi^ + [1(1 + 1) - 4] Y\^ = 0, (A5c) 



which are in a close correspondence with Eqs. (A3 1 
We begin with A^''\ which we decompose as 



(A6) 



in terms of harmonic components Am ■ There are 2Z -I- 1 
real terms in the sum, and the 2Z -I- 1 independent com- 
ponents of Ak-^k2---ki are in a one-to-one correspondence 
with the coefficients Am ■ Returning to the original rep- 



resentation of Eq. (|A2|) , we find after differentiation that 

DaA^'^ = m'XAak2 



■ki^'''^ ■ ■ ■ f^^^S and we conclude that 



rlrr 
'A 



(A7) 



This is the decomposition of Aa^ in vectorial, even- 
parity harmonics; this equation is valid for I ^ 
0. An additional differentiation using the last of 
Eqs. (|aT| reveals that DaDbA^'^ = -I^abA'^'^ + l{l 



I)ri^f)^yla6fe3 - fci ■ ■ ■ ■ From this we conclude that 



A 



(0 

AB 



\DaDi 



AB 



(A8) 



This is the decomposition of A^ab tensorial, even- 
parity harmonics; this equation is valid for / ^ {Oj !}• 

We next examine the odd-parity potentials. Their gen- 
eral structure is 



apql'i U dk3---k,l b 



and they satisfy the eigenvalue equations 
r^Y'^DcDdB^J^ + [l{l + 1) - 1]6W = 0, 



(A9a) 



(A9b) 



(AlOa) 
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0, 



(AlOb) 



the same as in the even-parity case. The potentials be- 
come 



B' 



(0 



A'^apq 



(Alia) 



■ n''' (Aiib) 



after transformation to angular coordinates. We wish 
to decompose these in the odd-parity harmonics of 
Eqs. (2.211, which satisfy the eigenvalue equations 



This is the decomposition of B^^ in vectorial, odd-parity 
harmonics; this equation is valid for I ^ Q. After a sec- 



ond differentiation we get 
l{l ^ m\eapqnPB\^^^ 
metrization of the indices we obtain 



and after sym- 



B' 



(I) 

AB 



l{l 

2 



1) 



{efDB + e^DA)DcB^'^ 



l{l - 1) 



\ ^ K7(i) -vim 



^AB- 



(A15) 



[i{i + i)~i\x'r 

AB-r[l{l + l)~^]X 



Im 
AB 



0, (A12a) 'pjjjg jg |-jjg decomposition of B^J^g in tensorial, odd-parity 
= 0, (A12b) harmonics; this equation is valid for I ^ {0, 1}. 



the same as in the even-parity case. 

We once more begin with B^^'^ :— Bk^...ki^^^ 
its decomposition 



■ fi'"' and 



(A13) 



Differentiating the first expression, multiplying 
this by the Levi-Ci vita tensor, and involving 



the second of Eqs. (All returns e^DsB'^'^^ 
elude that 



From this we con- 



B' 



'Db)B^'^ 



1 

7 ^ 



(A14) 



APPENDIX B: DETERMINANT OF THE 
HORIZON METRIC 



The horizon metric of Eq. (4.3 1 can be expressed as 
lAB — ^M^Q.AB +PAB- The metric determinant is given 



by ^ = 4M^ sin 6'(1 + - 



c 5 

gives 



1 A B 
A 



■ ) , where 



^^I^'^PCB and e — e\. Evaluating this 



^ = 4M^ sin 
1 



l + —M^{V"'+Qr 
45 ^ 



^AB\ 



0(5) 



(Bl) 



where the indices on ^^^d B\b are raised with Q'^^ . 
This convention will be used consistently below: all in- 
dices on tidal potentials will be lowered with ^Iab and 
raised with fl^^ . 
The identities 

' + 4.0", 
48 



cq Cqafc 
ab 


cq Cq AB 

^ '^AB'^ " 




cq K?q ab 
'^ab'^ 


cq Rq AB 


16 


gq gqafc 
ab 


Rq Rq AB _ 


-> 



7 



2" + 22" 



(B2a) 
(B2b) 
(B2c) 



can be established by direct computation, by making use 
of the definitions of the tidal potentials provided in Ta- 
bles |T] [n] [ml [TV] |V] and |VT] Inserting them into our ex- 
pression for ^7 returns the simple expression displayed 
in Eq. (|4.4|). 



APPENDIX C: 



CALCULATION OF TIDAL 
HEATING 



In this appendix we provide calculational details re- 
garding the heating of the black hole by the tidal in- 
teraction. The results derived here are used in various 
places in Sec. |IVD| 

To a degree of accuracy that is sufficient for our pur- 



poses, the inverse to the horizon metric of Eq. (4.3 I is 
1 



' 4 



6 



(f'^^+6^^^)+Af-2o(3). (CI) 



It is understood that on the right-hand side of this equa- 
tion, upper-case Latin indices are raised with f2'^^. 

We use "f^^ to raise indices on the shear tensor of 
Eq. (Ol. This yields 
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' 60 ^ 



E 



'oAB\ 



420 



45 
1 

36 



45 
1 

42 



2r,AB nA 



21' 



9 



63 



13 
126 



M-^0{6). 



(C2) 



To arrive at this result we made use of the identity 
Aac-^^b — ^^abAcdA^^ satisfied by any symmetric- 
tracefree tensor Aab on the unit two-sphere, as well as 
Eqs. (B2 ) from Appendbc|B We emphasize that while in- 
dices on aAB are raised with the physical horizon metric 
■y^^ , indices on all tidal potentials are raised with . 

To eval uate t he ri ght-hand side of Eq. (4.121 we use 
Eqs. (4.8 1 and (C2| to construct cr^so"*^, which we 



integrate over the horizon with the help of Eq. (4.4 1. 
The shear tensor is expressed as a multipole expansion, 
and its square consists of products of multipole mo- 
ments. Some of these products involve moments of the 
same order. For example, ctab'^^^ contains the term 
lM^{£\g + B\g){E'^^ + 3""^^), which is a product 
of quadrupole moments; such terms survive an angu- 
lar integration and contribute to the right-hand side of 
Eq. (4.121. Other products involve moments of different 



orders, and those integrate to zero. 

To evaluate an angular integral such as / £\b^'^ "^^i 
where dfi :— sin 9 dOdcj), we recall the definition £\b ~ 
^Ib^A^B and deduce the identity S^b^""^^ = Sl^S""''. 
We next substitute the expression for found in Table m 
and carry out the integration. These steps are repeatea 
for all other relevant products of multipole moments, and 
we obtain 



4tt j 


1 tj^B^^"^^ dn = 


^ c cab 

5 


(C3a) 


4nJ 




iBabB""', 

5 


(C3b) 


4tt , 


j t\j^t^^ dn^ 


Q 

c cabc 
^^^abc^ : 


(C3c) 


4ttJ 


f b^bB"^'^ dn = 


j^B^bcB . 


(C3d) 



All other integrat ions vanish, or contribute to the right- 
hand-side of Eq. (4.121 at order (M/TZ)^ and beyond. 



The general solution to Eq. (4.13) is 



Kq 

Stt 



A{v) = e««M(0) - Stt / J^(z;')e'"« dv' . (C4) 



After three integration by parts this becomes 



A{v) 



Stt 



1 



AiO)-T{0) T{0) 



+ T{v) + —T{v) + ^T{v) 

Kq Ki 



1 r d^ 

d^ 



^0 -'0 



(C5) 



The last term can be neglected, because K.Q^d^T/dv^ is 
of order {M /TVf and beyond the accuracy maintained in 
the computation of the flux function. The first collection 
of terms, those that depend on the initial conditions at 
w = 0, grow exponentially over a short time scale of or- 
der Kq ^ = 4M. We do not expect the black- hole area to 
behave in this way; we expect instead that the tidal in- 
teraction will produce a slow growth over a much longer 
time scale. To eliminate the run-away solution we de- 
mand that the solution satisfy the initial condition 



Stt 



1 



1 



A{Q)=T{Q) + + ^T{Q) + 0{Q). (C6) 



Under this restriction, Eq. (C5l reduces to 
1 



A{V) 



T{v) + —T{v) + ^T{v) + 0(9), (C7) 



and this is just Eq. (4.151. Notice that Eq. (C7l is com 



patible with the initial conditions of Eq. ( C6 1 . 

It is unusual, when dealing with horizons, to impose 
miiiaZ conditions on solutions to differential equations, as 



we have done in Eq. (C6 1. The reason, of course, is that 
event horizons are always identified by final conditions. 
We chose to proceed in this way because our horizon 
is not necessarily an event horizon, as we explained in 
Sec. |lV A] We also explained that the horizon becomes an 
event horizon when the tidal interaction switches off in 
the remote future. Under this condition we may impose 
the final condition that ^ = when v = oo. The exact 



solution to Eq. (4.13) is then 



A{v) = Stt / T{v')e-''°'^''' ""'^ dv' , 



(C8) 



and this leads once more to Eq. ( C7 ) after three integra- 
tions by parts. 



APPENDIX D: ALTERNATIVE DERIVATION OF 
THE BACKGROUND METRIC 

In this appendix we give a brief sketch of an alternative 
derivation of the background metric of Eqs. ( |3.3[ ). The 
derivation described in Sec. [V] relied on Zhang's observa- 
tion jl8j that the metric of a vacuum region of spacetime 
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around a timelike geodesic 7 is a functional of two sets 
of tidal moments £aia2---ai and Baia2 - ai', these are STF 
tensors that depend on proper time on the world line, 
and they are related to components of the Weyl tensor 
(and its derivatives) evaluated on the world line. Here we 
provide a precise definition of the light-cone coordinates 
and construct the metric systematic ally through order 
we recover the metric of Eqs. (3.3 1 and therefore 



confirm the validity of Zhang's observation. Our develop- 
ment follows the general methods developed in Ref. f2| . 



1. Formal definition of light-cone coordinates 



As in Sec. |II A| we consider a timelike geodesic 7 de- 
scribed by the parametric relations z°'{t), with r denot- 
ing proper time. On 7 we install an orthonormal tetrad 
(m" , e" ) of parallel-transported vectors. (We use primed 
indices to refer to points on the world line; unprimcd in- 
dices will refer to points off the world line.) For the time 
being we do not assume that 7 is situated in a vacuum 
region of spacetime; this assumption will be incorporated 
at a later stage. To assign light-cone coordinates {v,x°') 
to a point x off the world line we locate the unique future- 
directed null geodesic segment /3 that begins at x and 
ends at a point x' on the world line. (The construction 
requires that x be in the normal convex neighborhood 
of x'; the coordinates are defined in this neighborhood 
only. ) The advanced-time coordinate v is the value of the 
proper-time parameter r at this point: x' = z{t — v). 
And the quasi-Cartesian coordinates are defined by 



■.^-el,a^\x,x'), 



(Dl) 



where ej^, := 5 ga'p'S-'l and ct" := V" a is Synge's 
world function cr(x, x') [301 EI] differentiated with respect 
to its second argument. The points x and x' are related 
by the condition a{x,x') = 0, which indicates that the 
points are finked by a null geodesic segment. 
We define 



r := —an'U 



(D2) 



and state without proof that r is an affine parameter 
on (3; it decreases as the null geodesic approaches the 
world line. (The proof of this statement is contained in 
Ref. [2 .) The completeness relation — — + 

(5°^e" and the identity cr" ctq/ — 2a ^ imply that 
= Sabx""!^- It is useful to introduce 



x7r, (D3) 
and we use the completeness relation to write 

a"' = r-(u"' -f^"e^'). (D4) 



This is a decomposition of the displacement vector 
a" {x, x') in the tetrad (u" , e" ). The vector 

■■= <Ja/r (D5) 
is future-directed and tangent to null geodesic segment 
f3. Here cTq '■— VqCt is the world function (t{x, x') differ- 
entiated with respect to its first argument. 

Suppose now that the point x is moved to a neighbor- 
ing point X + 5x. The coordinates of the new point will 
he V + 5v and a;° + Sx"" , and to calculate the coordinate 
displacements we must locate the new point x' + 5x' on 
the world line, which is liked to x-f 5a; by a new geodesic 
segment (3 + 5(3. Using 5x°' = 5v and expanding 
a(x + 5x, x' + 5x') = to first order in the displacements, 
it is easy to show that 5v — —1^ Sx", so that 



(D6) 



The definition of the coordinates a;" in terms of the world 
function then implies that Sx"" — —{e'^,a"p,u^)6v — 

e°,cr"^ Sx^, so that 

d„x- - (eS,a>7.^') £^ - e''p,aC (D7) 



Here cr"^, := V" V^'cr is the second covariant derivative 
of the world function with respect to x' , while a^^ := 
Vq'CT denotes a mixed derivative with respect to each 
argument. 



2. Metric 

We begin with a computation of the inverse metric, 
with components 



vv 


g"^davdi3V, 


va 


g^^do.x'^dpv 


,ab 


g^^d^x'^dpx 



(D8a) 
(D8b) 
(D8c) 



Using Eq. (D6I and the fact that is a nu ll ve ctor we 
immediately find that g"" = 0. With Eq. ^% we get 



gva _ ^a^^^a^^ii^ simplify this with the help of 

Eq. JdsI and the identity ct^Ict'' = {a^af^) = V"'cr. 



the components 



With Eqs. (Dl) and (D3l, this is g"" = 0.°-. To obtain 



e^,ep,g '^cr ^cr ^ - (e„,cr ^,7/'^ jS2 



(D9) 



requires a much longer computation, and the result will 
be expressed as an expansion in powers of r. 
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We rely on the known expansions (see, for example, Ref. |72j ) 



1 ' A" 1 ' S' e' 



_ 0' 
<7a'/3 — 9 p 



1 ' 5' 1 ' S' e' 



7 



(DIO) 



(Dll) 



for the derivatives of the world function, in which g^p{x,x') is the parallel propagator [30l[7T]. We make the substi- 
tutions in Eq. (D9l and express cr" as in Eq. (D4|. 



The end result of a lengthy computation is 



5™ 


= 0, 




















gab 




o 




- p'-f^") - 




- P'^n'' 4 


- p^n'') - 




- g°r2^ ^ 






+ 




f p^r^'' - 






f Q°J7^ - 




+ 1.^(3.-- 


f s'^r^'' - 






+ 








45 ^ 






+ 0(r5). 







(D12a) 
(D12b) 



The inverse metric is expressed in terms of the potentials 

Pab = RaabO — {RacbO + RbcaCl)^'' + Racbd^'^^'^ — Pba 
Pa — RaOcQ^'^ ~ RacdQ^'^^'^ = Pab^^ , 
P = RcOdO^'^'' = Pan\ 
Qab = ~RaObO\c^'' + {RacbO\d + RbcaOld)^'^^''^ — Racbd\e^'^^'^^'^ — Qba 
Qa — ~ RaQcO\d^'^^'^ + RacdO\e^'^^'^^'^ — Qab^'' , 

Sab = RaObOlcd^'^^'^ ~ {RacbO\de + RbcaO\de)^''^'^^'^ + Racbd\ef^'^^'^^'^^'^ = Sba, 
Sa — RaOcO\de^'^^'^^'^ — RacdO\ef^'^^'^^'^^^ — Sab^^ , 
S — RcOdOlef^'^^'^^'^^'^ — Sa^"^, 

Uab — {RaOmoR"bcO + RbQmQR"acQ + RamboR'^cO + RbmaoR™co)^'^ 

\ f O t? E? I E? l?^^ I T? T?'^'- ODE? E> E?^ E? E?^^ 

+ {^J^aOcO^bOdO + J^amcO^ dbO ~r ^fcmcO^ daO ~ ^^oObO-ftOcOd ~ -n-amfcO^ cdO ~ ^bmaO^ cdO 

F? F? F?^ F? R'^^ F? F?'^ ^ O^O*^ 

^acmO-'L bdO ^bcmO^ adO ^acbrnJ^ QdO J^bcam^ OdO/ " ^' 

+ { — RaOcoRbdeQ — RbOcoRadeO + RacmdR"bea + RbcmdR^^aeO + RacboRdOeO 
+ RbcaClRdOeO + RacbmR"deQ + -Rbcam-R™eo) ^^^''^"^ — Uba^ 
Ua = {RaOmoR"cdO + RcmaoR"odo)^''^'^ ~ {RacmoR"deO + RacdmR"oeo)^'^^'^^'^ ~ Uab^^ , 
Vab = PacP% - PaPb = Ha, 
Va - PacP" - Pa/', 
V ^PcP" -P"^, 

I 



(D12c) 



(D13a) 
(D13b) 
(D13c) 
(D13d) 
(D13e) 
(D13f) 
(D13g) 
(D13h) 
(D13i) 



(D13j) 
(D13k) 
(D131) 
(D13m) 
(D13n) 



which are defined in terms of the frame components in Sec. |II A| and an overdot indicates differentiation 
of the Riemann tensor (and its derivatives) evaluated with respect to proper time. For example, Pab '■= 
on the world hne. We adopt the notation introduced RaObo — {Racbo + R-bcao)^'^ + R-acbd^'^^'^, where, for ex- 
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ample, Raobo '■= Ra' ^i' v' -x' e-a u^^i, u'^ ■ Notice that 
the derivative operator acts on the Riemann tensor only. 
The fact that the tetrad is parallel transported on the 
world line implies that the right-hand side can also be 
written as {Ra' n' p' u' u'^ );A'M^ , and we find that 
RaObo — dRaobo/dT- To avoid ambiguities with second 



derivatives of the Riemann tensor, we always differen- 
tiate with respect to proper time in the last step; for 
example, RaObo\c ■= ^QV/3V';7'A'e"'M'''e^ u^'ej'u^' . 



niques introduced in Sec. V A We obtain 



The inverse of Eqs. (D12) is calculated using the tech- 



Si; 



9va = ^a+ la 

1 

-15' 



-r^P, 



1 

12' 
1 



12 



Qc - i^r^P, - —r^Q 



15 



15 



gab = lab + la^lb 



1 , 1 



r^P, 



10 



r Qcd 



20 



r'^Scd 



— r'^Ucd+^r^Vcd + Oir^) 
20 45 I V ; 



r 



(D14a) 



(D14b) 



(D14c) 



where Tofc := Sab — ^a^b- We recognize the structure of 
Eqs. (5.7 1, with a clear decomposition of the metric into 



longitudinal and transverse pieces. 



Saci^bdSpq — SbqSdp) — Sad{SbcSpq ~ S^qScp) + 5aq{5bc5pd — 
Sbd^cp)- 

The frame components of the first spatial derivatives 
of the Weyl tensor are 



3. Decomposition of the Weyl tensor 

At this stage we demand that the Ricci tensor, its 
first derivatives, and its second derivatives, all vanish on 
the world line 7. This implies that the Riemann tensor 
Ra'p'-y'S' and its derivatives are equal to the Weyl ten- 
sor Ca'pi-y'S' and its derivatives. The symmetries of the 
Weyl tensor, the Bianchi identities, and the Ricci iden- 
tities then imply that the Weyl tensor (and its deriva- 
tives) can be expressed in terms of the tidal moments 
^ab: ^abc^ ^abcd^ Bab, Babe, and Babcd\ these were defined 
in Sec. 

The frame components of the Weyl tensor on the world 
line are given by 



CaObO — ^ab, 

CabcO — ^abpB^^, 

(~i cpq 

^ abed — ^abp^cdq^ 



(D15a) 
(D15b) 
(D15c) 



The last equation can also be written as 

Gabcd = Sac^bd — ^ad^bc — ^bc^ad + Sbd^ac, (D16) 

by making use of the general identity Sabp^cdq = 



CaObO\c — £ab\c, 
CabcO\d = ^abpB ^|^, 

Cabcd\e ~ ^abp^cdq^ 



(D17a) 
(D17b) 

(D17c) 



where £ab\c := CaQbo\c and Bab\c := h^'Pi^^loic These 
are related to the tidal moments by 



£ab\c = £abc + 2 i^acpB^b + '^bcpB^a) , (D18a) 
Bab\c = ^Babc - I {eacp£\ + efecp^'L) ■ (D18b) 

The frame components of the (symmetrized) second 
spatial derivatives of the Weyl tensor are 



CaObOKcd) — £ab\{cd), 
CnbrMidrA — ^abpB^ 



abcO\(de) — tabpO c|(tie)' 



c. 



abcd\{ef) — ~^abp^cdq^ 



/)' 



(D19a) 
(D19b) 
(D19c) 



ab\{cd) '■— C'oOf)OI(cd) and (cd) i^ap<?^^60| (cd) ' 



where £, 

These are related to the tidal moments by 



ab\{cd) 



bed 



B'bc + ^bcpB^d + ^bdpB^ac) 
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4 .. 1 .. .. .. .. 11 

+ ^^ab^cd — j{Sac£bd + ^ad^bc + ^bc^ad + Sbd^ac) + ^^cd^ab 

22 19. 20 

SabG,, + - [S^cGm + 5adGbc + 5,,Gad + kdGa.) + -5cdGab 



21 

4 

21 



21 



{^ab&cd + Sac^bd + ^ad&bc G, 



10. 



Bab\(cd) — 2 Babcd ^{^acp^' ijd 

1 



^bcp£^ad ^bdp£^ac 



^Sab^cd — J {SacBbd + SadBbc + hc^ad + ^bd^a 



11 



ScdBab 



— 2 {Sab^cd 



22 

21 

4 

21 



Babied) 
19 , 



'-.{SacBi 



bd 



SabHcd + ^ {SacHbd + SadHl 



be 



{SabScd + SacSbd + SadSbc)H. 



Bac£bd) 
\- 5bcHad 



21 

[SadBbc + BadSbc) 
20 



3 

SbdHa 



21 



ScdHab 



(D20a) 



(D20b) 



We introduced the notation 



Gab 


- £ £P ~ 


G 


— fPQf 


Hab 




H 


= 2£PWpg 



i~'api~' b' 

- BP''B, 



P9' 
CP 
ap^ by 



Metric in irreducible form 



(D21a) 
(D21b) 
(D21c) 
(D21d) 



K^{£ab£cd + BabBcdp^^n'n'n" 



£PB\n''n''n'=. 



+ Ccpgi-- a'-' b 



(D22) 



To write this in terms of irreducible tidal potentials we 
write a product of unit radial vectors such as riTi^fi^ri'^ 
as the STF decomposition 



The decomposition of the Weyl tensor and its deriva- 
tives in terms of the tidal moments is substituted within 

Pab^ Qab: Qa, Q, Sab, Sa, S, Uab, Ua, Vab, Va, and 

V. All of this is next substituted within our previous 
expression for the metric tensor. The manipulations re- 
quired to simplify the expressions are extremely lengthy 
(on the order of 35 pages of small handwriting) . We find 
that the terms that are linear in the Weyl tensor and 
its derivatives group themselves automatically into the 
irreducible tidal potentials introduced in Tables |T] and |TT] 
The organization of the terms that are quadratic in the 
Weyl tensor requires more work. For example, our initial 
expression for the quadratic terms in g^y is 

^quadratic ^ Lr^£p,£P^ - Bp.BP'l) 



21 



bBapB\)n''n'' 



= ^(afccd> _,_ l^ja6j^(cd> ^ac^{bd) gad^ibc) 
+ §bd^(ac) ^ gcd^{ab)\ 



15 



(D23) 



and we make the substitution within the metric function. 
Looking at the term SabScd^"'^''^'^^''' , for example, we 
obtain SabScd^^"'""^^ + ^SpaS^^'^"'''^ + ^£pq£P'', which can 
be expressed in the equivalent form S^ab^cd)^"'^''^''^'^ + 
^Spt^aS^'b)^"'^'' + Tb^pq^^'^- According to the definitions 
listed in Table 



III 



this is T"' + + f^V-", a superpo- 



sition of hcxadecapole, quadrupole, and monopole tidal 
potentials. 

Proceeding in a similar way with g^a and gab, we even- 
tually arrive at Eqs. (3.3 1. 
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